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ABSTRACT 

This  is  a  final  study  of  a  radiational  model  for  use  in  the  FNWC 
prediction  system.   This  model  utilizes  original  FNWC  data  over  oceanic 
gridpoints  which  were  adapted  to  the  FNWC  five-layer  initial  data 
analysis  in  four  mid-seasonal  studies  for  the  data  year  16  January 
1974,  16  April  1974,  16  July  1974  and  16  October  1973. 

The  radiational  model  parameterizes  the  large-scale  cloud  amounts 
in  two  layers.   The  primary  objective  was  to  tune  model  albedo  values 
to  those  taken  from  satellite  climatology  by  changing  the  cloud-reflec- 
tance coefficients  used  in  previous  studies.   The  albedo  values  were 
tuned  for  least  squares  deviation  relative  to  satellite  climatology 
albedos.   The  present  model  using  tuned  albedos  better  verified  the 
radiational  balance  at  the  top  of  the  atmosphere  against  satellite 
climatology  than  did  pre-existing  untuned  models,  both  on  a  seasonal 
and  annual  basis. 
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I.   INTRODUCTION 

This  thesis  is  a  refinement  of  previous  radiational  models  described 
by  (A)  Spaeth  (1975),  (B)  Meyers  (1975),  (C)  Beahan  (1975)  and  (D) 
Warner  (1974)  for  use  in  the  Fleet  Numerical  Weather  Central  (FNWC)  pre- 
diction system.   This  study  has  as  a  primary  objective,  the  comprehensive 
re-examination  of  the  radiational  physics  in  layers  comprising  the  ocean- 
atmosphere  system.   FNWC  atmospheric  soundings  defined  at  constant  pres- 
sure levels  for  the  four  mid-seasonal  dates  of  the  "data  year"  1973-74 
as  previously  examined  in  studies  (A,B,C,D)  were  utilized  in  this  study. 

The  application  of  the  FNWC  radiative  model  may  be  made  to  any  scale 
of  analysis  for  which  there  is  adequate  resolution  of  the  temperature 
and  moisture  data  in  the  vertical.   In  the  horizontal,  the  reliability 
of  the  data  used  here  is  consistent  with  that  of  the  FNWC  interpolation 
to  gridpoints  in  the  analysis  procedure.   Temperature  and  dew-point  data 
in  radiative  soundings  are  typically  reported  to  the  nearest  tenth  of  a 
degree.   The  radiative  computations  made  here  are  applied  to  FNWC  grid- 
points  and  are  designed  to  make  a  one-hour  forward-time  step  at  gridpoints 
in  the  FNWC  primitive  equation  forecast  model,  with  special  adaptions  to 
a-levels. 

The  specification  of  cloud  amounts  in  two  designated  layers,  one  at 
a  mid-level  and  the  other  at  a  low-level,  has  an  important  influence  on 
the  radiative-model  dispositions,  both  in  the  short-  and  long -wave 
spectral  regions.   Initially  in  (A,B,C,D),  the  specification  of  the 
fractional  amounts  of  CL  and  CL  had  been  based  on  large-scale  formulc 
tions  developed  by  Smagorinsky  (1960) ,  but  during  the  course  of  these 
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studies  it  was  found  more  realistic  to  modify  initial  CL-values  to 
CL'  =  2/3  CL.   This  reduction  in  cloud  amount  was  also  used  here  as  it 
had  been  chosen  to  prevent  CL '  from  exceeding  unity  and  to  afford 
better  agreement  with  satellite  climatology,  such  as  planetary  albedo, 
for  compatible  data  periods . 

With  the  reduced  cloud  coverages  CL  '  and  CL  '  from  the  earlier 
studies  (A,B,C,D),  it  was  possible  to  obtain  reasonably  close  agreement 
in  the  computed  terrestrial  net  flux  at  the  top  of  the  atmosphere  and 
that  observed  for  comparable  NIMBUS  III  subsatellite  points  and  data 
periods  (Raschke,  Von  der  Haar,  Bandeen  and  Pasternak,  1973).   However 
even  with  the  reduced  cloud  coverages  CL  '  and  CL  '  the  computed 
planetary  albedo  remained  generally  excessive,  particularly  in  tropical 
latitudes.   Hence  it  became  a  major  objective  of  this  particular  study 
to  modify  empirically  the  reflective  capability  of  the  cloud  layers. 
This  was  done  by  defining  a  general  factor  f  so  that  the  initial  choices 
(after  Rodgers,  1967)  of  cloud  reflectances  R  were  modified  to  R'  where 

R'  =  fR  o 

Systematic  substitution  of  the  cloud  reflectances  R,  wherever  they 
entered  the  solar  disposition  equations,  by  R1  then  led  to  a  relation- 
ship between  the  global  albedo  and  f.   Utilization  of  the  least  squares 
technique  to  minimize  the  differences  between  satellite  and  model  albedos 
over  a   geographic  sample  of  points  led  to  best-fit  value  of  the  "tuning- 
factor"  f.   Separate  values  of  f  were  deduced  by  least  squares  for  each 
season,  and  subselections  were  deduced  for  the  tropical  and  extratropical 
areas,  respectively. 
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The  modified  solar  cloud-reflectances  improved  the  agreement 
between  satellite  and  radiative  model  albedos  in  both  geographic 
areas  insofar  as  net  incoming  insolation  was  concerned.   The  terres- 
trial net  flux  at  the  top  of  the  atmosphere  was  unaffected  by  the 
choice  of  f ,  while  the  use  of  CL  '  and  CL  '  as  specified  gave  good 
agreement  with  satellite  terrestrial  net  flux  data  over  the  geo- 
graphic range  and  the  time-scales  concerned. 
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II.   DATA  PREPARATION 

A.   DATA  FIELDS 

1.   General  Considerations 

The  initial  temperature  and  humidity  data  used  in  this  study 
were  arranged  in  the  form  of  soundings  taken  along  four  oceanic  meri- 
dians (Fig.  1)  of  the  Fleet  Numerical  Weather  Central  (FNWC)  Northern 
Hemisphere  mid-seasonal  analyses  for  16  October  1973,  16  January  1974, 
16  April  1974  and  16  July  1974.   Oceanic  locations  for  these  computa- 
tions were  chosen  because: 

P 

(a)  constant  a-surfaces  (where  0   =  — )  of  the  FNWC  primitive 

equation  system  are  nearly  identical   constant  pressure  levels. 

(b)  The  maritime-area  soundings  are  more  likely  to  be  system- 
atically representative  of  the  set  of  zonally-distributed  gridpoints 
than  over  land. 

The  three  meridians  selected  over  the  Pacific  Ocean  were  lo- 
cated at  125W  (25  soundings),  170W  (25  soundings)  and  145E  (17  sound- 
ings) .   The  Atlantic  Ocean  meridian  was  located  at  35W  (26  soundings) . 
This  method  of  selecting  "soundings"  along  the  indicated  meridians 
of  the  FNWC  polar  stereographic  map  made  it  unnecessary  to  employ  spa- 
tial interpolation  between  original  data  gridpoints  along  the  meri- 
dians.  Data  along  line  3  in  the  Pacific  was  not  extended  southward  of 
gridpoint  (9,55)  because  they  fell  over  land  masses  (New  Guinea  and 
Northern  Australia)  where  the  surface  temperatures  and  other  sounding 
features  were  unrepresentative  of  the  oceanic  values. 
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Figure  1.   FNWC  polar  stereographic  grid  and  meridians  (lines  1,  2, 
3,  and  4)  selected  for  study.   The  longitudes  A  are  shown 
for  each  meridian  as  well  as  the  extent  considered  of  each 
'  meridian. 
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2.   Data  Treatment 

a.   Original  Soundings  and  Modifications 

The  gridpoint  soundings  were  taken  from  the  original  FNWC 
63-by-63  Northern  Hemisphere  analyses  of  T(p)  and  of  T  -  T   (the  dew 
point  depression)  at  standard  pressure  levels  up  to  and  including 
p  ~  100  mb.   Examples  of  such  original  soundings  were  shown  as  Table 
1(a)  in  Meyers  (1975,  p.  24).   Subsequently  each  original  FNWC  sound- 
ing was  transformed  — in  previous  studies  (A,  B,  C,  D)  of  this  series 
— into  what  has  been  termed  the  radiative  sounding  having  the  format 
shown  in  Table  I.   The  data  levels  of  the  radiative  sounding  contains 
essentially  the  five  FNWC  predictive  a-levels  (dotted  levels  in  Fig.  2) 

At  each  gridpoint  selected,  the  original  FNWC  humidity 

soundings  were  given  in  the  form  of  five  dew  point  depressions  over 

the  analyses  levels  from  925  mb  to  400  mb.   At  the  surface  (level 

k  =  10)  the  standard  instrument  level  vapor  pressure,  e  .  ,  was  trans- 

c        c  air' 

formed  into  the  surface  mixing  ratio,  q-,n/  by  means  of 

qin  =  621.97  (e  .  /1000)  .  (2-1) 

-LU  air 

To  obtain  radiative  soundings  as  in  Table  I,  it  is  neces- 
sary to  have  water  vapor  and  CO   absorber  masses  at  certain  required 
k-level  boundaries  (Fig0  2) .   All  radiative  soundings  in  this  study 
start  at  sea  level  with  the  approximation  of  surface  pressure 
it  =  1000  mb.   Therefore,  the  eleven  k-levels  correspond  closely  to  the 
FNWC  levels  P   =  1000.,  900.,  800.,  ....,  200.,  100.,  0.0  mb  and  in 
turn  toff,  =  1.0,  0.9,  ,  0.1,  0.0. 
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TABLE  I.   Example  of  a  radiative  sounding  at  gridpoint  (1,1)  for 
16  April  1974  with  mixing  ratio  listed  at  odd  k-levels  (Fig.  2) . 
Additionally,  water-vapor  and  C0_  absorber  masses  are  also  listed 
as  these  parameters  have  been  modeled  in  the  radiative  theory  pre- 
sented in  this  study. 


Mixing 

Pressure 

Temp 

Ratio 

(mb) 

(°C) 

(g/kg) 

Absorber  Masses 

Water  Vapor     CO? 

o  2 

(gm/cm^)        (cm/cm  ) 


1000 


25.60 


17.10 


900 

16.09 

12.08 

800 

11.21 

2.26 

45.53 

700 

5.00 

6.21 

600 

-2.10 

3.23 

83.53 

500 

-11.50 

2.57 

400 

-23.20 

1.35 

3.55 

113.35 

300 

-38.00 

0.65 

200 

-56.60 

3.61 

133.99 

100 

-80.70 

0.03 

3.61 

138.67 

0 

-80.70 

3.61 

143.35 
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Figure  2.   Five-layer  radiative  sounding  used  in  this  study.   Levels 
are  identified  by  their  values  on  the  k-scale,  while  layers  are  identi- 
fied by  their  level  boundary  indices  in  parentheses,  e.g.  (8,10). 
Pressure-scaled  water  vapor  and  CO  mass  increments  Au  and  AC,  respec- 
tively are  integrated  relative  to  the  surface  and  the  resulting  U  and  C 
are  carried  at  even  levels.   The  temperature  T  is  retained  at  all  levels, 
Amounts  of  clouds  CLj   and  CL2   in  the  layers  shown  have  been  para- 
meterized for  consideration  of  their  radiative  effects. 
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b.  Radiative  Temperature  Profiles 

The  gridpoint  temperatures  were  listed  at  each  mandatory 
level  of  Table  I  between  1000.,...,  100  mb  (i.e.,  between  k  =  10,..., 
1) .   The  temperature  was  assumed  to  be  isothermal  from  100  mb  to 
0.0  mb.   The  temperature  at  level  k  =  10  was  set  equal  to  the  FNWC 
listed  sea-surface  temperature.   The  radiative  sounding  temperatures 
for  the  remaining  k-levels  were  obtained  from  either  their  correspond- 
ing listed  temperature-level  or  by  a  three-point  Lagrangian  interpola- 
tion scheme  [Eq.  (2-1) ,  Spaeth,  1975] ,  to  level  k  when  the  listed 
FNWC  temperature  profile  did  not  include  the  value  T  . 

c.  Radiative  Moisture  Profiles 

Similarly,  the  moisture  profiles  of  Table  I  have  been  ob- 
tained by  an  interpolative  procedure  over  the  original-level  FNWC  mix- 
ing ratios  to  those  required  at  k-levels  in  a  manner  analogous  to  that 
discussed  by  Spaeth  (1975,  pp.  29-31). 

d.  Pressure-Scaled  Absorber  Masses 

The  pressure-scaled  water  vapor  and  the  carbon-dioxide 
scaled  absorber  masses  were  calculated  for  the  six  even  numbered  k- 
levels  (Fig.  2)  using  Eqs.  (2-8,  2-9,  2-10  and  2-11)  respectively  as 
outlined  by  Spaeth  (1975,  pp.  31-33) .   These  equations  use  essentially 
the  mixing  ratios  of  water  vapor  and  of  CO_  at  odd  k-levels. 

B.   CLOUD  PARAMETERIZATION 

The  relative  humidities  (RH)  at  levels  k  =  5  and  k  =  9  are  used 
in  the  calculations  of  the  fractional  cloud  covers  CL  and  CL   in 
layers  (4,6)  and  (8,9)  respectively  (Fig.  2).   The  equations  for 
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parameterization  of  the  two  fractional  cloud  amounts  are  as  follows: 

CL   =  2/3  [2.0  (RH(5))  -  0.7]  (2-2a) 

CL2  =  2/3  [3.33(RH(9))  -  2.0]  (2-2b) 

The  bracketed  part  of  the  equations  (after  Smagorinsky,  1960)  were 
reduced  by  the  1 /3  factor  in  an  attempt  to  tune  cloud  amounts  to  obtain 
albedo  values  in  closer  agreement  with  the  recent  satellite  radiational 
climatology  of  Raschke  et  al.,  (1973).   The  "2/3  CL"  parameterization 
set  forth  in  Eq.  (2-2)  is  used  here  for  estimating  large  scale  radia- 
tional effects  only.   Thus  small-scale  convective  activity  and  a   priori 
climatological  effects  were  not  considered  in  specifying  the  form  of 
Eq.  (2-2) .   Tuning  of  the  model-albedo  values  of  this  study  by  varying 
cloud  reflectance  coefficients  with  respect  to  season  and  latitude  will 
be  discussed  in  Section  IV. 

C.   CLOUD-AREA  COVERAGES 

Fractional  cloud  amounts,  CL  and  CL  ,  were  computed  at  each  grid- 
point  for  levels  k  =  5  and  k  =  9  by  Eqs.  (2-2a,b) .   In  addition  the 
gridpoint  area  may  be  thought  of  as  broken  into  random  fractional  seg- 
ments of  size 

W(0,0)  =  (1-CL1)  d-CL2)  (2-3a) 

wherein  there  is  a  combination  of  clear-over-clear  segments  in  the 
layers.   Similarly,  the  gridpoint  area  has  the  fractional  area  of  cloud 
coverage 


W(l,l)  =  CL   *CL2  (2-3b) 


2  2 


of  an  upper  overcast  amount  overlying  a  lower  overcast.   Likewise  the 
area-combinations  of  overcast  over  clear,  and  clear  over  overcast 
areas  in  two  layers,  may  be  visualized  as  occurring  with  the  weights 


and 


W(1,0)  =  CL   *(1-CL2)  (2-3c) 


W(0,1)  =  (1-CL  )  *CL2  (2-3d) 


respectively . 

For  radiational  computations  it  was  useful  to  carry  the  relative 
weights  or  fractions  of  the  gridpoint  area  exposed  to  the  specified 
cloud-layer  combinations.   Henceforth,  the  symbols  denoted  by  W(0,0), 
W(l,l),  W(1,0)  or  W(0,1)  indicate  the  fractionally  overcast  (1)  or 
clear  (0)  cloud-area  combinations  in  the  indicated  layers  (Fig.  2), 
with  the  first  index  1  or  0  referring  to  layer  CL  ,  k  =  5,  and  the 
second  to  CL  ,  k  =  9. 

The  usefulness  of  the  cloud-area  weighting  device  will  be  clarified 
in  Sections  III  and  IV,  where  the  procedures  for  the  terrestrial  and 
solar  radiational  computations  are  discussed  and  the  results  are 
summarized  over  the  set  of  soundings. 

A  measure  of  the  effective  cloud-cover  area  which  has  been  found 
useful  in  previous  radiational  studies  has  been  the  total  opaque  cloud 
cover,  CL,  referring  to  the  amount  of  thick  cloud  cover  overhead  re- 
gardless of  the  level.   For  the  cloud  model  presented  here  CL  may  be 
expressed  as 


CL  =  CL   +  CL   -  CL   *CL   .  (2-4) 
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III.   TERRESTRIAL  RADIATION 

A.   THEORETICAL  AND  EMPIRICAL  BASIS 

Empirical  formulas  were  developed  by  Sasainori  (1968)  for  flux 
emissivities  in  the  atmosphere  associated  with  computations  for  the 
radiative  balance  requirements  of  the  NCAR  General  Circulation  Model. 
Sasamori  derived  the  empirical  emissivity  formulas  for  water  vapor 
and  CO  by  comparison  with  the  theoretical  values  built  into  the 
Yamamoto  Radiation  Chart  (1952) .   The  Yamamoto  chart  has  proved  to 
be  quite  accurate  for  numerical  checks  of  the  Sasamori  emissivities. 
This  chart  was  also  used  in  the  previous  studies  (A,  B,  C,  D)  as  a 
systematic  guide  for  integration  of  the  radiative  transfer  formulas 
developed  by  Martin  (1972,  1975),  who  adapted  the  Sasamori  emissiv- 
ity formulas  to  the  particular  layers  of  interest  in  the  gridpoint 
computations  of  the  FNWC  primitive  equation  model  (Fig.  2). 

The  essential  long-wave  (IR)  net-flux  parameters  required  for 

use  in  this  study  are  the  following: 

* 
F    =  IR  net  flux  at  earth,  k  =  10 

* 

F0   =  IR  net  flux  at  level   k  =  8 
o 

* 

F.    =  IR  net  flux  at  level  k  =  6 

6 
* 

F„    =  IR  net  flux  at  level   k  =  4 
4 

* 

F    =  IR  net  flux  at  level   k  =  2 

In  addition  the  IR  net-flux  divergence  coolings  to  be  computed  at 

each  gridpoint  are 

F810  =  IR  net-flux  divergence  in  the  layer  (8,10) 
F68   =  IR  net-flux  divergence  in  the  layer  (6,8) 
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F46  =  IR  net-flux  divergence  in  the  layer  (4,6) 
F24   =  IR  net-flux  divergence  in  the  layer  (2,4). 
In  the  Radiation  Balance  Studies  in  the  series  A,  B,  C,  D  only 
F610  and  F26  were  computed  because  time  restraints  in  the  present 
FNWC  operational  heating  package  have  prevented  the  use  of  greater 
resolution  in  the  vertical.   Here,  four  flux  divergences  are  computed 
in  order  to  examine  more  closely  the  variability  of  the  flux  diver- 
gences over  the  layer  thicknesses  reduced  to  approximately  200  mb 

each.   To  compute  these  four  flux  divergences  it  was  necessary  to 

*       * 

utilize  additional  formulas  for  F   and  F    as  developed  by  Martin 

o         4 

(1975) . 

In  order  to  make  IR  net-flux  calculations  along  the  path  of  in- 
tegration, there  must  be  a  physically  sound  representation  of  the 

emissivity  (£   )  as  a  function  of  both  water  vapor  and  CO^  absorber 
wc  2 

masses  in  layers  along  the  sounding.   For  a  complete  discussion  of 
the  emissivity  formulas  used  in  the  quadrature  scheme,  refer  to 
Spaeth's  Appendix  A  (1975). 

B.   NET  FLUX  FORMULATIONS 

The  radiative  sounding   as  depicted  in  Table  I  was  computed  as 
the  combination  of  parameters  U(k,10),  C(k,10)  and  T   for  each  re- 

.K 

quired  level,  k  =  10,  8,  ...  1,0.   Cloud  parameters  CL   and  CL  were 
also  computed  by  Eq.  (2-2)  at  each  gridpoint  and  in  general  are  both 
non-zero.   The  grid  area  was  then  considered  to  be  composed  of  areal 
fractions  (weights)  defined  in  Eqs.  (2-3a,b,c,d)  and  denoted  by  the 
symbols  W(0,0),  W(l,l),  W(1,0),  W(0,1). 

The  composite  net  flux  F     (CL  ,  CL  )  at  level  k  =  10  at  each 
gridpoint  is  then  constructed  by  using  the  appropriate  weight  factors 
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to  multiply  the  reference  net  flux  F    computations  defined  for  the 
four  special  cloud-cover  cases  already  defined  in  Section  II. C: 

Fio*(0'0)'  Fio*(1'0)'  Fio*(0'1)'  ^o**1'11 

It  therefore  follows  that 

F10*(CL1'  CL2)  =  W(0,0)F10*(0,0)  +  W(1,0)F10*(1,0) 

+  W(0,1)F1()  (0,1)  +  W(1,1)F10*(1,1)  .      (3-1) 

* 
The  reference  net  fluxes  F    of  Eq.  (3-1)  are  associated  with  (1) 

clear  skies  in  both  layers,  (2)  overcast  in  the  upper  layer  only, 

(3)  overcast  in  the  lower  layer  only  and  (4)  overcast  in  both  layers 

respectively. 

Spaeth  (1975)  has  listed  these  reference  net  flux  formulations 

in  his  Eqs.  (3-6),  (3-7)  and  (3-8).   Using  the  definitions  of  W(0,0), 

W(1,0),  W(0,1)  and  W(l,l),  F   *  (CL  ,  CL  )  can  be  shown  to  assume  the 

form 

F10*(CL1'  CL2)  =  [1"CL2  H(B10-B6)-.5[ewc(8r10)(Bl£)-B8) 
+  (ewc(8,10)  +  Bwc(6/10))(B8-B6)]} 

+  (1-CLo   )  (1-CL-,   ){B-.5[(e   (6,10) 

*         '6       wc 

+  ewc(4,10))(B6-B4)  +  (Cwc(4,10) 
+  Ewc(2/10))(B4-B2)  +  (£wc(2,10) 

+  £wc(i,o))(B2-Bl)  +  ewc((o,io),T1)*B1]} 

+  CL9  {(B.  -B_)  [1.5G   (9,10)]}  .  (3-2) 

^     10   9      wc 


Here 


B.  =  1.170403  x  10  ?  T,  4  (3-3) 

k  k 
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is  the  Stefan-Boltzmann  blackbody  flux  in  langlies  per  day. 

Further,  £   (U,  ,C,  ,10)  is  the  combined  water-vapor  and  CO,,  emissiv- 
wc   k   k  2 

ity  along  the  path  from  level  10  to  level  k.   This  emissivity  is  con- 
sidered by  Sasamori  to  be  temperature  independent  for  T  >   210K,  whereas 

£    represents  the  temperature  dependent  emissivity  applicable  for 
wc 

T  <  210K.   [See  pp.  136-137,  Spaeth  (1975);  Sasamori  (1968)]. 

* 
The  formulas  for  F   (CL  ,  CL  )  ,  k  =  2,4,6,  and  8  have  been  developed 


by  Martin  (1975)  in  a  manner  analogous  to  the  derivation  of  the  weighted 
1 
10 


* 

F    .   The  results  are  reproduced  as  the  following  equations 


FQ   =  [1-CLJ  {b-.5[£   (6,8) (B-B,) 
8  1     8      wc        8   6 

+  (£   (6,8)  +  £   (4,8))  (B-B.) 
wc  wc         6   4 

+  (£   (4,8)  +  £   (2,8))  (B  -BJ 
wc         wc        4   2 

+  (G   (2,8)  +  £   (1,8)) (B  -B.)  +  £   ((0,8),T1)*B.]} 

wc         wc        2      1     v/c        1    1 

+  CL   [1-.5E   (6,8)]  (BQ-B.)+  CLn  (1-CL_)  [1- .  5£   (8, 10)]  (B.  -BQ) 
1        wc         8   6      1       2        wc         10   8 

+  (1-CL  )  (1-CL.)  [l-.5£   (8,10)]  (Bin-BQ)  (3-4) 

J.         Z  WC  ±U   o 

F/  =  [1-CL^  {B8-.5[£wc(6,8)(B8-B6)  +  ^  (4,6)  (Bg-B^ 

+  (£wc(4'6)  +  ewc(2'6))(B4"B2)  +  (£wc(2'6) 
+  Ewc(lf6))(B2-B1)  +  £wc((0,6),T1)*B1]} 

+  (1-CL_)  (1-CL_){(B1 n-B_)  [l.-.5(£   (6,8) 

X         /.  1U   O  wc 

+  £   (6,10))]}  +  CLn  {(B0-B.)* 
wc  1     8   6 

[1-.5E   (6,8)]}  +  CL.  (1-CL.)  {(Bln-BQ)* 

wc  1      2      10   8 

[l-.5(£   (6,8)  +  £   (6,10))]}  .  (3-5) 

wc         wc 
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F4*  =  [1-CL1]  {B8-.5[ewc(4,6)(B6-B4) 

+  (Ewc(4'6)  +  £wc(4'8))(B8-V  +  EWC(2'4)(B4-B2) 

+  (Ewc(2'4)  +  Ewc(1'4,)(VBl)  +  5wc((0.4)fT1)*B1]} 

+  CL   {B4-.5[£wc(2,4) (B4-B2)  (3-6) 

+  (£   (2,4)  +  £       (1,4))(B  -Bn)  +  £       ((0,4)fT.)*B.]} 

wc         wc   '      2   1     wc        11 

+  (1-CL. )  (1-CL_)  {(B.  -B_)  [l-.5(£   <4,8)  +  £       (4,10))]} 
1      1  10   o        wc         wc 

F2*  =  [1-CI^]  {B8-.5[£wc(2,4)(B4-B2)  +  (^(2,4) 
+  Ewc(2f6))(B6-B4)  +  (£wc(2,6)  +  £wc(2,8))* 

(VV  +  £wc(1'2)(B2-Bl)  +  5wc((Of2)#T1)*B1]}         (3-7) 

+  (1-CL.)  (1-CL.)  {(B.  -B_)  [l-.5(£   (2,8) 
1       2      10   o         wc 

+€       (2,10))]}  +  CL.   {B-.5[£   (2,4)(B-B_) 
wc  1     4      wc        4  2. 

+  £       (1,2MB  -B.)  +  £       ((0f2)fT1)*B.]}  . 
wc       2   1     wc        11 

As  was  described  by  Spaeth  (1975,  Section  III.B.5.)  concerning  the 

use  of  the  composite  case,  Eq.  (3-2)  can  be  reduced  to  give  expressions 

* 
for  F    for  the  various  reference  cloud-cover  cases  (0,0),  (1,1),  (1,0) 

*       * 

and  (0,1).   The  resulting  schematics  in  the  case  of  F   and  F   are 

o        4 

depicted  as  the  unhatched  area  in  Figs.  3(a,b,c,d)  and  4(a,b,c,d), 

*     *       * 

respectively  below.   Similar  graphs  for  F    ,  F    and  F   can  be  found 

on  pp.  41-43  of  Spaeth  (1975) ,  and  are  not  reproduced  here. 

* 

A  typical  gridpoint  listing  of  the  IR  net-flux  computations,  F 

(k  =  10, 8,... 2),  has  been  reproduced  in  Table  II  for  gridpoint  (1,1) 
based  upon  the  radiative  sounding  of  16  April  1974  (see  Table  I) .   The 
printout  procedure  involves  computation  of  the  reference  net-flux 
values 

Fk*(0,0),  Fk*(l,0),  Fk*(0,l),  Fk*(l,l) 
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from  each  of  Eqs.  (3-2),  (3-4),  (3-5),  (3-6)  and  (3-7)  for  level 
k  =  10,  8,... 2  respectively.   Then  Eq.  (3-1)  with  the  appropriate 

weight-factors  of  Eqs.  (2-3a,b,c,d)  has  been  utilized  to  derive  the 

* 

composite  F   (CL  ,  CL  )  values  that  are  listed  on  the  bottom  line  of 

Table  II. 


TABLE  II.   A  sample  listed  of  IR  net-flux  computations,  weighting 
factors  and  composite  values,  F   (CL  ,  CL  ) ,  as  computed  for  grid- 
point  (1,1)  for  16  April  1974.   Net  flux  values  in  ly  min-1. 

Cloud  Case   Weight 


(CL   CL2) 

W(CL  ,  CL2) 

Fio 

F8 

F6 

F 

4 

F2 

(0,0) 

.1002 

.1664 

.2288 

.2649 

.3160 

.3487 

(1,0) 

.1497 

.0943 

.1331 

.0915 

.1706 

.2338 

(0,1) 

.3007 

.0512 

.1563 

.2341 

.2871 

.3202 

(1,1) 

.4494 

.0512 

.0586 

.0607 

.1706 

.2338 

* 

P,  -compc 

site  values 

.0692 

.1171 

.1379 

.2202 

.2713 

C.   TROPOSPHERIC  COOLING  BY  LAYERS;  F     COOLING 

For  each  mid-seasonal  day   listed  and  at  each  gridpoint,  an  IR  net- 
flux  computation  in  the  format  of  Table  II  is  easily  converted  into 
four  sets  of  layer  cooling  effects. 


*  * 

F810  =  F    -  F 

8      10 

*  * 

F68   =  F.   -  F0 
6      8 

*  * 

F46   =  F.   -  F„ 

4      6 

*  * 

F24   =  F„   -  F„ 
2      4 


-1 


(3-8) 


These  layer  cooling  rates  (ly  min  "  have  then  been  collected  in  meri- 
dional cross-section  format  for  each  longitude  under  study  (see  Section 
V)  and  by  season. 
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The  overall  tropospheric  cooling  rate  by  IR  net  flux  is  then  given 

*      * 
simply  by  F   -  F    at  each  gridpoint  for  the  date  of  the  radiational 

sounding.   The  tropospheric  cooling  rates  computed  are  then  identical 

to 

*      * 
F    -  F     =  F24  +  F46  +  F68  +  F810  .         (3-9) 

*       * 

The  values  of  F   -  F    so  deduced  are  discussed  on  both  a  seasonal 

and  a  zona lly-ave raged  basis  in  Section  V. 

It  will  suffice  to  discuss  here  the  zonally-averaged  values  of 
F     (CL  ,  CL„)  as  computed  by  the  long-wave  radiational  model  previously 
presented  in  this  section,,   These  results,  listed  simply  as  F    in 

Table  III,  will  be  discussed  as  a  function  of  seasons,  latitudes  and 

* 

CL  (total  opaque  cloud  cover  given  by  Eq.  (2-4)).   The  listings  of  F 

in  Table  III  are  essentially  as  extracted  from  computations  in  the 

* 

format  of  Table  II  followed  by  meridional-averaging  of  F    across 

constant  latitude  lines.   Finally  the  zonally-averaged  annual  values 

* 
of  both  F    and  CL  have  been  computed  by  arithmetic-averaging  over 

the  four  mid-seasonal  results  at  each  five-degree  increment  of  latitude 

from  20S  to  65N  (cf.,  Eq.  (5-1)). 

* 

The  model-annual  values  of  F    presented  in  Table  III  are  presented 

with  those  derived  from  Budyko  (1956) ,  which  in  turn  are  listed  in  the 
final  column  of  Table  III.   Corresponding  values  of  total  opaque  cloud 

cover,  CL,  for  the  Budyko  climatology  were  not  available  so  that  only 

* 

a  general  comparison  of  the  two  annual  F    zonally-averaged  distribu- 
tions is  possible. 
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Table  III  depicts  the  zonally-distributed  values  of  F    and  of  CL 

at  the  earth's  surface.   In  a  seasonal  comparison  of  the  model-computed 

*  * 

F    values  it  is  clearly  shown  that  F    (CL)  is_  a_  decreasing  function 

of  cloud  cover.   There  is  a  clear-cut  tendency  in  each  season  for  a 

* 
maximum  value  of  F    to  be  located  in  the  subtropics  (latitudes  15N- 

* 

25N)  .   Also  there  is  evidence  of  a  high  latitude  (55N-60N)  minimum  F 

associated  with  a  concentration  of  maximum  cloud  cover  CL.   An  outstand- 
ing variation  is  the  transition  in  the  Southern  Hemisphere  latitudes 
(20S  -  10S) ,  which  has  small  cloud  cover  in  local  winter  and  compari- 
tively  large  cloud  cover  during  the  other  three  data  periods.   This 
cloud-cover  variation  corresponds  in  general  to  the  ITCZ  behavior  in 
these  latitudes  across  the  indicated  seasons;  so  that,  F    is  a  maximum 
in  mid-July  and  a  relative  minimum  in  the  period  January-April. 

The  final  entry  in  Table  III,  namely  "Wt.  Avg.",  is  the  Northern 

Hemisphere  mean,  cosine-weighted  with  respect  to  latitude  (cf.,  Eq.  (6-1)) 

* 

The  cosine-weighted  F    values  show  a  minimum  in  Northern  Hemisphere 

summer  with  no  clear-cut  differences  in  CL.   This  summer  minimum  is 
attributable  to  higher  downward  flux  with  increased  summer  vapor 
pressures. 

D.   OUTGOING  IR  NET  FLUX  TO  SPACE 

1.   Parameterization  Formula  for  Top-of-Atmosphere  IR  Net  Flux 

In  the  four  previous  radiation  studies  (A,  B,  C,  D) ,  an  approxi- 
mation to  the  IR  net  flux  to  space,  designated  as  FF2,  was  computed 
from  the  radiative  sounding  at  each  gridpoint.   FF2  was  essentially 
an  extrapolation  of  F    to  the  top  of  the  atmosphere  obtained  by  delet- 
ing the  downward  IR  flux  due  to  stratospheric  water-vapor  and  CO  .   A 
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more  precise  expression  for  the  net  IR  flux  to  space  was  introduced  here 
after  Martin  (1975) c   The  development  is  analogous  to  the  quadrature 
summations 


B=0 

£   dB  (3-10) 

wc 


'-■  f 


Bio 

for  the  various  reference-cloud  combinations  (0,0),  (1,0),  (0,1)  and  (1,1) 
of  Fig.   5(a,b,c,d)  respectively.   The  final  quadrature-formula  is  then 
obtained  as  the  weighted  net-flux  result  as  was  also  done  in  Eqs .  (3-4) , 
(3-5),  (3-6),  (3-7)  and  is  listed  below. 

FT  -  [l-CL1]{B8-.5[(£wc(0,l)  +  Cwc(0,2))(B2-B1) 

+  (£wc(°'2)  +  £wc(°'4))(VB2) 

+  (e   (0,4)  +  £   (0,6))  (B-B.)  +  (£   (0,6)  +  E   (0, 8)  )  (B-B,.)  ]  } 
wc         wc        6   4      wc         wc        o   6 

+  CLn  {B-.5[(£   (0,1)  +  £   (0,2)) (B-B.)  (3-11) 

14      wc         wc        2   1 

+  (£   (0,2)  +  £   (0,4))  (B  -BJ]} 
wc         wc        4   2 

+  (1-CLj  (l-CL.){[l-.5(£   (0,8)  +  £   (0,10))]  [B   -B_]} 

1        Z  WC  WC  1U   o 


where  CL  ,  CL  are  given  by  Eq.  (2-2).   Note  that  F   of  Eq.  (3-11)  is 
representable  by  the  unhatched  areas  of  Figs.  5(a,b,c,d)  for  the  various 
reference-cloud  cases  and  that  F  has  no  downward  IR  flux  corresponding 
to  the  level  k  -  0. 

2.   Comparisons  of  F  with  Mid-Seasonal  Satellite  Climatology 

Comparison  was  made  of  computed-model  values  of  F  with  satellite 
measurements  of  total  long-wave  flux  to  space  (after  Raschke  et  al.,  1973) 
for  the  Nimbus  III  mid-seasonal  periods  most  nearly  comparable  to  that 
of  the  FNWC  data.   F   for  each  mid-seasonal  date  was  computed  over  the 
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four  meridians  considered.  These  F  values  were  then  averaged  across 
the  four  meridians  to  get  a  mean  zonal  distribution  of  the  type  shown 
in  Table  IV. 

Table  IV  shows  the  zonally-averaged  model  values  compared  with 
those  extracted  from  Raschke.   Raschke's  results  were  obtained  by 
averaging  across  the  same  oceanic  meridians  as  those  used  in  this  study. 
Again,  in  the  bottom  line  of  each  column  in  Table  IV  is  listed  the 
cosine-weighted  mean  of  each  set  of  column  values  for  the  Northern 
Hemisphere  only. 

The  zonally-averaged  values  computed  by  the  F  -model  are  very 
close  to  those  reported  by  Raschke,  especially  between  0-65N  and  in 
the  Northern  Hemispheric  means.   The  limitations  of  the  comparisons 
made  here  are  obvious,  when  it  is  recalled  that  between  latitudes  20S- 
5S  and  between  60N-65N  there  are  fewer  than  four  meridional  lines 
available  for  computing  the  listed  zonal  values  in  Table  IV.   For  all 
other  zonally-averaged  values,  four  meridional  lines  were  used  in  the 
averaging. 

The  close  comparison  between  model-values  of  F   from  (3-11)  and  those 
essentially  derived  from  satellite  climatology  tend  to  support  the 
cloud  parameterization,  Eq.  (2-2) ,  insofar  as  IR  net  flux  is  concerned. 
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IV.   SOLAR  RADIATION 


A.   COMPOSITION  OF  SOLAR  INSOLATION 

At  the  top  of  the  atmosphere  (k=0)  this  study  assumed  a  solar  con- 
stant of  2.00  ly  min  "  (Joseph,  1971).   Furthermore,  this  constant  was 
assumed  subject  to  a  four  percent  attenuation  above  the  tropopause  due 
to  ozone  and  oxygen.   Thus  the  effective  solar  constant  at  level  k=2 
in  this  study  is  1.92  ly  min 

To  compute  the  effective  solar  insolation  at  the  tropopause  the 
following  formula  was  used 

F(2)  =  S[— ]~2  Cos  Z  (4-1) 

m 


where   S  is  the  effective  solar  constant  at  level  k=2  and 

r/r  =  ratio  of  the  actual  earth-sun  distance  to  the  mean 
m 

earth-sun  distance,  a  function  of  the  Julian  date 
Cos  Z-   cosine  of  the  zenith  angle,  a  function  of  the 
Julian  date  determined  by 

Cos  Z  =  Sin  <$>   Sin  6  +  Cos  <j)  Cos  6  Cos  h  .  (4-2) 

The  symbols  on  the  right  side  of  Eq.  (4-2)  are  defined  as  follows: 
$  =  latitude 
6  =  solar  declination 
h  =  hour  angle 
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The  Smithsonian  Meteorological  Tables  (List,  1958)  list  the  ratio  of 

r/r  and  the  solar  declination,  6,  for  the  mid-seasonal  dates  appli- 
m 

cable  to  this  study  and  reproduced  in  Table  V. 


TABLE  V.   Values  of  the  ratio  of  the  earth-sun  radius  vector,  r/r  , 
and  of  the  solar  declination  angle,  5,  used  in  this  study  for  year 


1973  - 

1974. 

Date 

r/r 
m 

6 

16  January 

0.98372 

21.07917°S 

16  April 

1.00333 

8.48333°N 

16  July 

1.01644 

21.50000°N 

16  October 

0.99717 

8.22500°S 

The  value  of  sin  9  was  calculated  using  one  of  two  different  formulas, 
depending  on  the  data-line  used  for  the  computations,  in  terms  of  the 
FNWC  map  coordinates  (I, J)  as  in  Eq.  (4-3a,b) .   Conversely  for  these 
lines  one  may  solve  for  I  in  terms  of  sin  (p   as  in  (4-3c,d)  : 


■>    o  a        o •   a    973.752  -  2 (32-1 ) 2  , .    _  , 

Lines  1,3,4   Sin  9  =  ■ —  (4-3a) 

973.752  +  2(32-1) 

t-            e-   a    973.752  -  (32-I)2  ,.  ... 

Line  2        Sin  <p  = (4-3b) 

973.752  +  (32-1) 

Lines  1,3,4    I  =  32  -  22.065  [-  Cos ,  ^   ]  (4-3c) 

1  +  Sm  q) 

Line  2        I  =  32  -  31.205  [-  Cos ,  ^   ]  (4- 3d) 

1  +  Sin  9 

'1,...,25      for  Lines  1,2 
I  =(8 ,...,25      for  Line  3 


63,.., 38      for  Line  4 
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Here  I  is  the  abscissa  distance  on  the  FNWC  grid  (Fig.  1)  and  varies 
by  line  as  described  in  Section  II.   The  soundings  for  lines  1,  2,  and 
3  were  all  taken  at  OOOOGMT  with  the  solar  noon  existing  at  the  180th 
meridian;  therefore  the  hour  angles  for  these  three  lines  were  55°, 
10°,  and  35°,  respectively.   For  line  4,  the  soundings  were  taken  12 
hours  earlier  with  solar  noon  at  the  Greenwich  meridian,  giving  an 
hour  angle  for  line  4  of  35°. 

A  very  simple  partition  of  solar  insolation  was  utilized  in  this 
study  after  Joseph  (1971) .   It  consisted  of  dividing  the  insolation 
F(2)  into  two  parts  at  level  k-2.   One  part  was  considered  to  include 
all  wavelengths  X  >  .9  ym  where  absorption  by  water  vapor  and  carbon 
dioxide  bands  are  the  most  prevalent  attenuation  processes  in  clear  air. 
This  part  of  the  solar  spectrum  was  termed  the  F(A)  energy  and  con- 
sidered subject  to  water-vapor  absorption  but  not  to  Rayleigh  scatter- 
ing.  For  shorter  wavelengths,  X  <    .9  ym,  absorption  of  the  solar  in- 
solation energy  by  water  vapor  was  considered  negligible.   This  part  of 
the  solar  insolation  was  denoted  F(S)  suggestive  of  the  fact  that  it 
was  subject  only  to  Rayleigh  scattering  attenuation  in  clear  air.   The 
two  partitions  are  formulated  after  Joseph  (1971)  as  follows: 

F(A)  =  .349  F(2)  (4-4) 

F(S)  =  .651  F(2)  .  (4-5) 

In  this  study,  the  introduction  of  two  cloud  decks  produced  cloud-reflec- 
tivity effects  upon  both  the  F(A)  and  F(S)  solar  energy  insolations. 
However,  in  the  clear  areas  around  any  gridpoint  only  the  absorption- 
attenuation  applies  to  the  F (A)  insolation,  while  only  Rayleigh  scatter- 
ing-attenuation applies  to  the  F(S)  insolation. 
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B.   DISPOSITION  OF  F(S)  INSOLATION 

In  the  disposition  of  the  F(S)  insolation,  Joseph  (1971)  determined 
that  Rayleigh  scattering  reflectance  to  space  by  clear  skies  (after 
Coulson,  1959)  could  be  effectively  approximated  by  least  squares  in 
the  following  form 

a(R)  =  .085  +  .25074  [log  (—-)  Sec  z)  ]  (4-6) 

o 


where  P   =  1013.25  mb.   In  Eq.  (4-6),  tt/P   =  1  in  view  of  the  fact  that 
o  o 

mean  sea  level  pressure  7T  is  close  to  1000  mb.   Also 


Sec  z  -  (Cos  z) 


with  Cos  z  given  by  Eq.  (4-2) . 

The  surface  albedo  a(G)  is  another  reflective  parameter  utilized  in 
this  study.   Over  oceanic  areas  the  following  formula  for  a(G)  after 
Gates  et  al  (1971) ,  was  utilized: 

a(G)  =  max  {.06,  .06  +  .54  (.7  -  Cos  z) }  .  (4-7) 

As  described  in  Section  III,  four  combinations  of  reference-cloud 
cases  are  possible  with  a  two-layer  cloud  model.   The  disposition  of 
F(S)  under  each  of  these  cases  will  be  discussed  in  the  remainder  of 
this  subsection. 

1.   Clear  Sky  Case 

In  the  clear  sky  (0,0)  case  the  F(S)  insolation  was  subjected 
to  both  Rayleigh  scattering  reflectance  a(R)  and  the  surface  reflectance 
a(G).   Considering  the  likelihood  of  a  succession  of  multiple  reflections 


42 


between  earth  and  atmosphere,  the  F(S)  insolation  actually  penetrating 
the  earth's  surface  after  scattering  is  given  by 

isio(o,o)  -  F(S)  [l-a(R)]  [l+a(R)a(G)  +  ...  (a(R)a(G)  )n 

+  ...]*(l-a(G))  (4-8a) 

that  is,  by 

IS10(0,0)  =  F(S)  [l-a(R)]  [l-a(G)]/[l-a(R)a(G)  ]  (4-8b) 

2.   Cloudy-Sky  Cases 

In  the  three  cases  in  which  clouds  were  present,  F(S)  insolation 
absorbed  by  the  ground  at  each  gridpoint  was  computed  using  the  follow- 
ing equation  (after  Arakawa,  1972) : 

IS10(1,1)  =  F(S)  (l-R(D)  (1-R(2))  (l-a(G)) 

*{l-[R(l)R(2)  +  R(2)a(G)"+  R(l)a(G) 

-  2R(l)R(2)a(G)]}"1  .  (4-9) 

As  indicated  by  the  notation  (1,1),  denoted  CL  =  CL   =  1.0,  Eq.  (4-9) 
is  the  formula  used  in  calculating  F(S)  insolation  absorbed  by  the 
ground  in  the  case  where  overcast  clouds  are  present  at  both  levels  of 
Fig.  2.   Also  in  Eq.  (4-9),  initial  values  of  cloud-reflectance  were 
chosen,  namely  R(l)  =  .54  for  the  mid-level  clouds  between  k=4  and  6, 
and  R(2)  =  .66  for  the  low-level  clouds  between  k=8  and  9.   Both  cloud- 
reflectance  values  are  as  suggested  by  C.  D.  Rodgers  (1967) . 

For  all  the  other  cloud  cases,  the  following  changes  were 
applied  to  Eq.  (4-9).   In  the  (1,0)  case  (CL  =  1.0,  CL  =  0.0),  the 
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desired  earth-absorbed  insolation  is  obtained  by  setting  R(2)  =  0.0 
in  (4-9) ,  from  which  it  follows  that 

IS10(1,0)  =  F(S)  (l-R(l))  (l-a(G))/[l-R(l)a(G)]  .  (4-10) 

In  the  case  (0,1),  one  sets  R(l)  =  0.0  in  (4-9)  so  that  (4-9)  simpli- 
fies to 

IS10(0,1)  =  F(S)  (1-R(2))  (l-a(G))/[l-R(2)a(G)].  (4-11) 

Note  that  with  a  cloud  overcast  present,  the  Rayleigh  clear-sky  scatter- 
ing a(R)  does  not  appear  in  Eqs .  (4-9),  (4-10)  or  (4-11),  but  is  included 
empirically  in  the  cloud  reflectances  R(l)  and  R(2) . 
3.   Composite  F(S)  Insolation  at  Earth 

Equations  (4-8) ,  (4-9) ,  (4-10)  and  (4-11)  were  utilized  in  the 
computation  of  the  cloud-weighted  F(S)  insolation  penetrating  the 
earth ' s  surface  considering  the  areal-weights  of  the  cloud  combinations 
denoted  by  (0,0),  (1,1),  (1,0)  and  (0,1)  about  a  gridpoint.   The  re- 
sultant F(S)  insolation  penetrating  the  earth's  surface  denoted  by 
IS10  is  therefore  expressible  as 

IS10(CL  ,  CL  )  =  IS10(0,0)  W(0,0) 
+  IS10(1,1)  W(l,l) 
+  IS10(1,0)  W(1,0) 
+IS10(0,1)  W(0,1)  .  (4-12) 

Here  the  weighting  factors  W(0,0),  W(l,l),  W(1,0)  and  W(0,1)  are  com- 
puted in  Eqs.  (2-3a,b,c,d)  respectively.   Note  finally  that  the  part 
of  F(S)  insolation  reflected  to  space  is  found  by  subtracting  IS10 
(CL  ,  CL  )  from  F(S)  . 


44 


Table  VI  lists  the  F (S) -disposition  resulting  from  a  particular  radia- 
tive sounding  at  gridpoint  (1,1)  on  16  April  1974.   The  individual 
computations  of  IS10  as  they  apply  for  the  possible  overcast-clear 
layer  cases  are  made  under  the  heading  "IS10."   The  difference 

REFS  =  F(S)  -  IS10  (4-13) 

in  each  case  represents  F (S) -insolation  reflected  to  space  while 

STRAN  =   :j:S1°  (4-14) 

l-a(G) 

has  been  computed  as  that  portion  of  the  F (S) -insolation  incident  at 
the  sea  surface  just  prior  to  transmission  by  the  surface.   Note  that 
no  absorption  in  air  has  been  included  in  the  computations  of  Table  VI, 
and  that  the  only  absorption  permitted  is  that  implicit  in  IS10. 
Finally  at  the  bottom  of  each  column,  e.g.,  IS10,  the  composite  value 
has  been  computed  by  means  of  the  weighting  scheme  of  Eq.  (4-12) . 

TABLE  VI.   A  sample  listing  of  values  of  F(S)  insolation  (ly  min   ) 
computed  at  gridpoint  (1,1)  for  16  April  1974  using  Eqs.  (4-6,..., 
4-14) . 


"loud- 

-case 

Weight 

(CI^, 

CL2) 

W(CLX, 

CL2) 

IS10 

REFS 

STRAN 

(0,0) 

.0993 

.4305 

ol751 

.5216 

(1,0) 

.1484 

.2539 

.3517 

.3076 

(0,1) 

.3016 

.1921 

.4135 

.2327 

(1,1) 

.4507 

.1340 

.4716 

.1696 

Compo: 

=ite- 

F(S) 

values 

.2014 

.4041 

.2441 
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C.   DISPOSITION  OF  F(A)  INSOLATION 

The  fractional  portion  of  the  solar  insolation  subject  to  absorption 
by  atmospheric  water-vapor  and  carbon  dioxide  is  covered  in  the  follow- 
ing subsections. 

1.   Clear-sky  Case  (0,0) 

The  Manabe-Moller  absorptivity  function  provided  the  necessary 
absorptivity  values  for  the  key  layers  in  this  case.  The  form  of  this 
absorptivity  function  is 

.303 
a(2,k)  -  .271[U(2,k)  Sec  z]      .  (4-15) 

Here  a  is  the  absorptivity  applied  to  the  pressure-scaled  water  vapor 
mass  between  levels  2  and  k  (Fig.  2)  along  the  zenith  slant-path  angle 
z.   The  resultant  absorbed  insolational  energy  in  the  particular  layer 
(2,4)  is  then  given  by  the  Manabe-Moller  relation 

303 
A24  -  0.271F(A) [U(2,4)  Sec  z] *     .  (4-16) 

In  the  same  manner  A26,  A28  and  A210  are  found.   Then  the  absorbed  in- 
solation in  the  layers  (4,6),  (6,8)  and  (8,10)  are  computed  by 

A46  =  A26  -  A24  (4-17) 

A68  =  A28  -  A26  (4-18) 

A810=  A210-  A28  .  (4-19) 

Water-vapor  mass  above  level  2  was  assumed  negligible  in  the  F(A)  dis- 
position of  the  solar  insolation. 

By  subtracting  A210  from  F(A),  the  direct  transmission  of  F(A) 
insolation  impinging  at  the  earth's  surface  was  determined.   The  trans- 
mission of  F(A)  insolation  is  then  further  reduced  by  the  transmissivity 
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(l-a(G))   after  surface-reflectance, which  leads  to  the  earth-absorbed 
insolation 

303 
IA10(0,0)  =  F(A){l-.271[U(2,10)Sec  z]  *    }(l-a(G)).      (4-20) 

The  transmitted  energy  impinging  upon  the  earth  just  prior  to  absorp- 
tion is 

TRANA(0,0)  =  IA10(0,0)/[1-CC(G)]  .  (4-21) 

In  the  remainder  of  this  subsection,  which  discusses  the  cloudy 
layer  cases,  representative  cloud  reflectivities  and  cloud  absorptivi- 
ties  were  initially  adopted,  after  C.  D.  Rodgers  (1967) ,  for  the  two 
possible  cloud  layers.   These  initial  cloud  reflectivities  were  RA(1)  = 
.46  and  RA(2)  =  .50  while  the  cloud  absorptivities  were  A(l)  =  .20  and 
A (2)  =  .30.   Note  that  the  reflectivities  for  the  F (A)  wavelengths  are 
somewhat  smaller  than  those  adopted  for  the  F(S)  wavelengths.   The  pro- 
cedure of  considering  the  cloud  conditions  to  be  overcast  whenever  they 
appear  and  then  applying  the  appropriate  weighting  factors  in  the  com- 
posite summation  will  again  be  followed  as  in  Sec.  IV. B. 

To  simplify  the  discussion  for  the  cloud-covered  cases,  the  (1,1) 
case  will  be  presented  first  as  it  contains  representative  type  equa- 
tions for  the  remaining  two  cases,  (1,0)  and  (0,1). 
2.   Overcast  in  Both  High-  and  Low-cloud  Layers 

The  following  set  of  formulas  illustrate  the  model   disposition 
of  incoming  solar  insolation  (F(A))  from  level  k=2  to  the  earth's  sur- 
face and  permits  determination  of  the  amount  of  insolation  absorbed  by 
the  atmospheric  layers  and  by  the  earth's  surface.   The  dashed  separa- 
tion lines  are  introduced  to  subdivide  the  absorption  and  reflection 
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physics  of  the  model  into  subsections  which  permit  the  analysis  to 
proceed  more  or  less  within  successive  200  mb  layers 0   The  equations 
relate  to  the  parameters  in  Fig.  6  where  the  insolations,  A24,  A46, 
A68,  A89  and  A910,  etc.,  represent  the  contributions  to  the  insola- 
tion absorbed  in  the  layers  involved.   Symbols  F2,  F4,  F6,  F8  and  F9, 
etc.,  depict  the  streams  of  insolation  passing  through  the  indicated 
level.   A  vertical  arrow  implies  the  direction  of  insolation  passage, 
i.e.,  y  denotes  downward  insolation,  t  upward-reflected  insolation, 
and  4-i  downward-reflected  insolation.   Terms  involving  the  symbols 
"TD",  as  expressed  by  the  functions  of  Eq.  (4-22),  indicate  the  Manabe- 
Moller  transmissivities  for  diffuse  insolation  in  the  layer  beneath  an 
existing  cloud.   In  the  latter  situation,  the  term  Sec  Z  in  Eq.  (4-15) 


is  effectively  replaced  by  the  mean  slant-path,  Sec  Z  —  5/3  (Katayama, 
1966),  e.g., 

303 
TD68   =  1-.271[U(6,8)  5/3] 

TD910  =  l-.271[U(9,10)5/3]  *3°3  .  (4-22) 

In  accordance  with  the  above  definitions  and  as  described  by  Fig.  6, 
the  formulas  for  the  (1,1)  case  are  listed: 

303 
F4I  =  F(A)  (1-.271[U(2,4)  Sec  z]     ) 

F4+  =  F44-(RA(1)) 

303 
F2+  =  F4+  (1-.271[U(2,4)  Sec  z]     ) 

A24  =  F(A)  -  F4y  +  F4+  -  F2t  (4-23a) 
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A46    =   F4+(A(1))  (4-23b) 


F64-    =  F44-    -   F4+    -   A46 

F8i    =  F6l(TD68) 

F8t    =  F84-(RA(2)) 

F6t    =  F8+ (TD68) 

F6^+=  F6t (RA(1) ) 

F8l4-=  F64-4-  (TD68) 

A68   =  F6+    -   F8+   +   F8+    -    F6+    +   F644    -   F8-H  (4-23c) 


A89   =    (F8+    +    F844)  (A(2)  ) 

F94-    =   F8+    -   F81    +   F84-4-    -   A89 

F10l=    F94-(TD910) 

FlOt=   FlOHa(G)) 

F9+    =   FlOt (TD910) 

F9H=   F9t(RA(2) ) 

F10H=   F944(TD910) 

A910   -   F94-    -   F104-   +   FlOt    -   F9t    +    F9H    -   F1044 

A810   =  A89   +   A910  (4-23d) 

Note  that  the  effect  of  multiple  reflections  between  clouds  or 
between  the  earth's  surface  and  the  lower  cloud  has  been  incorporated 
to  include  only  the  effect  of  two  reflections,  with  the  lowermost  re- 
flecting surface  absorbing  the  remaining  impinging  insolation.   Computa- 
tions indicated  that  the  insolation  remaining  after  two  reflections  was 
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F(A) 
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^F2t 


A24+ 


F64-4-  /      lF6i 
A68I4-    A68  + 


A68+ 


A244- 


A91044   A910+ 
/    F10I  \ 


/  F10I g(G) 


F10+- 


IA10 


Figure  6.   Schematic  representation  of  F(A)  insolation  disposition  in 
the  case  of  two  overcast  layers. 
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too  small  to  warrant  the  consideration  of  further  reflections.   Also 
insolation  reflected  upward  from  a  lower  interface  (cloud  or  ground) 
to  the  base  of  an  upper  cloud  deck  has  not  been  subjected  to  upper 
cloud  absorption.   This  tends  to  reduce  very  slightly  the  secondary 
cloud-absorption. 

From  Eq.  (4-23d) ,  the  impinging  F (A)  insolation  at  the  earth's 
surface  can  be  expressed 

TRANA(1,1)  =  F10+  +  F1044  (4-24) 

The  F(A)  insolation  which  is  actually  absorbed  by  the  earth's  surface 
(see  Fig.  6)  may  be  written  as 

IA10(1,1)  =  FlO!(l-a(G))  +  FlO^i  (4-25) 

3.   Disposition  of  F (A)  Insolation  with  an  Upper  Overcast  Only 

With  a  single  cloud  layer  present  only  at  the  upper  level  (Fig. 
7)  the  equations  depicting  the  model-disposition  of  incoming  insolation 
becomes  a  simplified  subset  of  the  previous  case: 

F44-  =  F(A)  (1-.271[U(2,4)  Sec  z]*3°3) 

F4+  =  F44-(RA(1)) 

F2+  =  F4t (1-.271[U(2,4)  Sec  z]  ) 

A24  =  F(A)  -  F4^  +  F4f-F2t  (4-26a) 

A46  =  F4+(A1))  (4-26b) 
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Figure  7.   Schematic  representation  at  F(A)  insolation  disposition 
with  an  upper  overcast  layer  only. 
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F64-  =   F44-    -   F4t    -   A46 

F8l  =   F64(TD68) 

F104-  =   F64-(TD610) 

FlOt  =   FlOl(a(G)  ) 

F8t  =   FlOt (TD810) 

F6t  =    F10+(TD610) 

F644  =   F6t  (RA(1)  ) 

F844  =   F644(TD68) 

FlO+4-  =   F8I4-(TD810) 

A68  =  F64-  -  F84-  +  F8t  -  F6t  +  F644  -  F8++  (4-26c) 

A810   =  F84-  -  F10+  +  FlOt  -  F8t  +  F844  -  F1044  (4-26d) 

The  variables  used  above  are  defined  in  a  similar  manner  to  those  in  the 
(1,1)  case.  The  impinging  insolation  at  the  earth's  surface  may  be  for- 
mulated as  follows 

TRANA(1,0)  =  F10+  +  F104-4-  (4-27) 

while  the  F(A)  insolation  absorbed  by  the  earth  in  this  case  is  given  by 

IA10(1,0)  =  FlO+(l-a(G))  +  F10+I  .  (4-28) 
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4.   Disposition  of  F (A)  Insolation  with  a  Low  Overcast  Only 

With  an  overcast  lower  cloud  layer  the  model-disposition  symbols 
are  as  depicted  in  Fig.  8  and  are  physically  related  as  follows: 

303 
F4+    =  F(A) (1-.271[U(2,4)  Sec  z]     ) 

.303 
F6i    =  F(A) (1-.271[U(2,6)  Sec  z]     ) 

303 
F8I    =  F(A)  (1-.271[U(2,8)  Sec  z] *    ) 

F8t    =  F8HRA(2)  ) 

F6+    =  F8t(l-.271[U(6,8)  Sec  z] *  °3) 

303 
F4  +    =  F8t(l-.271[U(4,8)  Sec  z] '    ) 

.  303 
F2t    =  F8t(l-.271[U(2,8)  Sec  z]     ) 

A24    =  F(A)  -  F44-  +  F4t  -  F2f  (4-29a) 


A46  =  F4I  -  F64  +  F6f  -  F4  +  (4-29b) 

A68  =  F6I  -  F84-  +  F8i  -  F6+  (4-29c) 

A89  =  F8I  (A(2)) 

F94-  =  F84-  -  F8+  -  A89 

F104-  =  F9l(TD910) 

FlOt  =  F10l(a(G) ) 

F9t  =  Fl0t(TD910) 

F9ii  =  F9t(RA(2)) 

F104-4-  =  F94-KTD910) 

A910  =  F9l  -  F104-  +  FlOi  -  F9+  +  F9-R  -  F1044 

A810  =  A89  +  A910  (4-29d) 
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Figure  8.   Schematic  representation  of  F(A)  insolation  disposition 
with  a  lower  overcast  layer  only. 
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Again  the  variables  are  defined  as  stated  before  in  the  (1,1) 
case.   Likewise  the  incident  flux  is  defined  at  the  earth's  surface  as 

TRANA(0,1)  =  Fl(H  +  FIOH  (4-30) 

while  that  portion  which  is  absorbed  by  the  earth  (Fig.  8)  is 

IA10(0,1)  =  FlO+(l-a(G))  +  F10I+  .  (4-31) 

In  Eqs.  (4-25),  (4-28),  and  (4-31),  the  quantity  F1044  is  small  enough 
in  each  case  so  that  no  further  reflections  from  the  earth  were  con- 
sidered. 

5.   Composite  F(A)  Layer-Absorptions  and  Surface-Absorption 
Insolation 

As  has  been  previously  discussed,  the  standard  grid-area  weight- 
ing scheme  of  this  study  was  applied  to  obtain  composite  values  of  the 
absorbed  F (A) -insolation  in  key  layers  and  also  within  the  earth's  sur- 
face.  The  weighting  factors  applied  to  the  corresponding  overcast- 
combination  absorption  quantities  provided  the  following  composite 
results: 


A24(CL  ,CL  )  =  A24(0,0)  W(0,0)  +  A24(l,0)  W(1,0) 

+  A24(0,l)  W(0,1)  +  A24(l,l)  W(l,l)         (4-32) 


A46(CL  ,CL2)  =  A46(0,0)  W(0,0)  +  A46(l,0)  W(1,0) 

+  A46(0,l)  W(0,1)  +  A46(l,l)  W(l,l)         (4-33) 


A46(CL  ,CL  )  =  A68(0,0)  W(0,0)  +  A68(l,0)  W(1,0) 

+  A68(0,l)  W(0,1)  +  A68(l,l)  W(l,l)         (4-34) 
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A810(CL  ,CL  )  -  A810(0/0)  W(0,0)  +  A810(l,0)  W(1,0) 
+  A810(0,l)  W(0,1)  +  A810)l,l)  W(l,l) 


(4-35) 


IA10(CL  ,CL  )  =  IA10(0,0)  W(0,0)  +  IA10(1,0)  W(1,0) 

+  IA10(0,1)  W(o,l)  +  IA10(1,1)  W(l,l)  . 


(4-36) 


The  weighting  factors  W  (0,0)  ,  . .  .  ,  (Wl,  1)  were  listed  in  Eqs  .  (2-3a,b,c,d) , 
and  A24(0,0),  A46(0,0),  A68(0,0),  A810(0,0)  and  IA10(0,0)  are  given  in 
each  clear  sky  case  (0,0)  about  each  gridpoint  by  Eqs.  (4-15) ,  (4-16) , (4-17) 
(4-18),  (4-19)  and  (4-20)  respectively. 

The  results  for  the  absorption  in  layers  (2,4),  (4,6),  (6,8), 
(8,10)  and  at  the  surface,  level  k  =  10,  are  shown  in  the  following  table: 


TABLE  VII.   A  sample  listed  of  F(A)  insolation  values  (ly  min   ) 
computed  at  gridpoint  (1,1)  for  the  16  April  case. 


Cloud-Case   Weight 


CL1'  CL2 

W(CL1#< 

(0,0) 

.0993 

(1,0) 

.1484 

(0,1) 

.3016 

(1,1) 

.4507 

Composite-F (A)  values 


A24  A46  A68  A810  IA10  REFA  TRANA 

.0454  .0356  .0386  .0418  .1348  .0285  .1634 

.0633  .0559  .0280  .0156  .0482  .1137  .0581 

.0464  .0381  .0728  .0771  .0239  .0664  .0286 

.0633  .0559  .0432  .0285  .0104  .1233  .0126 

.0564  .0485  .0494  .0425  .0324  .0953  .0391 


In  the  computational  scheme  indicated  by  the  entries  of  Table 
VII,  the  reflected  F(A)  insolation  to  space  has  been  depicted  by  the 
symbol  REFA,  and  its  values  follow  from 


REFA  =  F(A)  -  A24  -  A46  -  A68  -  A810  -  IA10 


(4-37) 
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whereas  the  TRANA  dispositions  are  given  by  Eqs.  (4-21),  (4-24),  (4-27) 
and  (4-30)  or  by  its  weighted-mean  value  in  the  case  of  TRANA-composite, 
6.   Absorptivity  (ABA)  by  Layers 

Here  the  (fractional)  absorptivity  as  well  as  the  actual  insola- 
tion values  absorbed  in  the  layers  are  considered.   In  the  computation 
of  absorptivity,  which  is  fractional  absorption,  the  total  undepleted 
insolation  at  the  top  (k=0)  is  used  as  a  base.   The  following  equation 
was  utilized  in  this  calculation: 

-2 
FADJ  =  2.00 (r/r  )    Cos  z  .  (4-38) 

m 

The  absorptivity  of  the  troposphere  ABA  was  computed  from  the  ratio  of 
the  insolation  absorbed  in  the  troposphere  to  the  insolation  incident 
at  the  top  of  the  atmosphere  rather  than  at  k=2: 


A24  +  A46  +  A68  +  A810  (A    ^,a^ 

aba  =  ^rr^; *  Jyj 

FADJ 


D.   ALBEDO  (ALB)  OF  THE  EARTH-TROPOSPHERE  SYSTEM 

In  considering  the  planetary  albedo,  the  reflected  insolations  of 
the  earth-troposphere  systems  in  both  the  F (A)  and  F(S)  insolational 
regions  must  be  recalled  by  the  program.   Thus  REF  is  computed  at  each 
gridpoint  as  the  sum  of  the  reflected  insolation  energy  in  F (A) ,  denoted 
REFA  in  Eq.  (4-37)  and  the  reflected  part  of  F(S)  denoted  REFS  in  Eq. 
(4-13)  : 

REF  =  REFS  +  REFA  (4-40) 

Finally  the  planetary  albedo  is  related  to  FADJ  through 

RPF1 
ALB  =.—  .  (4-41) 
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E.   COMPOSITE  ABSORPTIVITY  (ABG)  BY  THE  EARTH-SURFACE; 
COMPOSITE  ATMOSPHERIC  TRANSMISSIVITY  (ATRAN) 

1.  Absorptivity  (ABG)  of  Earth 

By  summing  the  weighted  values  of  F(S)  and  F(A)  portions  of 
the  incoming  insolation  entering  the  earth,  the  total  insolation  ab- 
sorbed at  the  earth's  surface  was  computed.   This  quantity  when 
divided  by  the  extraterrestrial  insolation  gave  the  fractional  absorp- 
tivity (ABG)  of  the  earth's  surface.   The  equation  for  ABG  was 

m„„    IA10  +  IS10  ,.  ... 

ABG  = (4-42) 

FADJ 

where  IA10,  IS10,  and  FADJ  were  defined  previously  by  Eqs .  (4-36), 
(4-12)  and  (4-38)  respectively. 

2.  Transmissivity  (ATRAN)  of  the  Troposphere 

Also  computed  was  the  total  insolational  energy  TRAN,  incident 
at  the  earth's  surface  just  before  absorption  by  the  surface.   This 
calculation  is  given  by 

TRAN  =  TRANA  +  STRAN  .  (4-43) 

Here  STRAN  =  [IS10/  (1-ot  (G)  ]  was  previously  defined  in  Eq.  (4-14)  and 
a(G)  was  given  in  Eq.  (4-7).   TRANA  has  also  been  defined  as  the  weighted 
value  of  TRANA  (0,0),  TRANA  (1,1),  TRANA  (1,0)  and  TRANA  (0,1)  given  by 
Eqs.  (4-21),  (4-24),  (4-27)  and  (4-30).   Note  also  in  justification  of 
STRAN  that  the  four  cases  for  IS10  of  (4-8) ,  (4-9) ,  (4-10) ,  and  (4-11) 
each  have  the  common  factor  (l-a(G))  in  the  numerator  and  therefore  each 
transmitted  F(S)  insolation  component  available  at  the  earth  just  before 
absorption  needs  only  be  divided  by  (l-a(G)).   TRAN  may  thus  be  viewed 
as  the  total  insolational  energy  incident  at  a  pyrheliometer  located  at 
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earth.   The  (fractional)  transmissivity  of  the  troposphere  (ATRAN)  is 
then  computed  from 

ATRAN  =  TRAN/FADJ  .  (4-44) 

Note  finally  that  the  major  dispositions  of  the  total  insolation  at 
the  indicated  map  times  have  now  been  identified  by  the  fractional 
values,  ALB,  ABA  or  ABG,  and  ATRAN,  representing  the  reflectivity 
(albedo),  absorptivity  of  air  or  earth,  or  atmospheric  transmissivity 
as  the  case  may  be. 

F.   ALBEDO  TUNING  BY  COMPARISONS  WITH  SATELLITE  CLIMATOLOGY 
1.   General  Remarks  Concerning  a  Need  for  Tuning  Albedos 

In  the  four  previous  radiational  studies  (A,  B,  C,  D)  compari- 
sons were  made  between  the  solar-insolation  albedo-model  (ALBMOD)  com- 
putations and  the  satellite-climatology  albedo  (ALBRAS)  of  Raschke 
et  al.  (1973).   The  results  of  these  comparisons  indicated  excessively 
high  values  of  ALBMOD  especially  in  the  tropical  and  subtropical  oceanic 
areas.   The  question  was  raised  by  the  previous  investigators  as  to 
whether  vertically-structured  convective  cloud  elements  in  the  tropics 
would  have  as  high  a  reflective  capability  as  attributed  to  the  large- 
scale  cloud  masses  existing  primarily  in  horizontal  layers  as  indicated 
by  the  parameterization  of  Eq.  (2-2) .   It  should  be  recalled  that  the 
initial  cloud-reflectances  used  were  as  suggested  by  C.  D.  Rodgers 
(1967)  and  are  listed  in  Sees.  IV. B. 2.  and  IV.C.l. 

For  the  purpose  of  tuning  the  model  solar  dispositions,  and 
in  particular  ALBMOD  with  ALBRAS,  the  initial  values  of  Rodgers'  cloud- 
reflectances  were  simply  adjusted  by  a  multiplicative  factor  (f)  which 
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turns  out  by  the  least-squares  fitting  to  be  smaller  than  unity.   This 
device  allows  a  greater  fraction  of  solar  radiation  to  penetrate  down- 
ward through  the  cloud- types  in  both  the  tropical  and  subtropical  areas, 
and  provides  better  agreement  between  ALBMOD  and  ALBRAS. 
2.   Method  of  Tuning  ALBMOD  to  ALBRAS 

The  model-tuning  process  involved  several  steps.  First,  the 
new  variable  f  was  defined  and  allowed  to  range  in  value  from  1.2  to 
0.25  by  0.05  intervals.  The  f-values  were  then  multiplied  by  the 
initially  used  cloud-reflectance  values  (after  Rodgers,  1967)  giving 
twenty  sets  of  new  cloud-reflectivity  values.  For  each  set  of  new 
values,  ALBMOD  and  subsequently  the  difference  A (ALB)  between  ALBMOD 
and  ALBRAS ;  i.e., 

A (ALB)  =  ALBMOD  -  ALBRAS  (4-45) 

was  recomputed  at  each  gridpoint  for  the  four  midseasonal  dates. 

Because  cloud-types  differ  in  structure  over  tropical  and  extra- 
tropical  oceanic  areas  the  resulting  A (ALB)  data,  from  (4-45),  were 
divided  into  two  populations  denoted  TR  and  HL.   The  population  TR  in- 
cluded the  set  of  gridpoints  located  south  of  and  including  25N,  while 
HL  included  the  set  of  gridpoints  located  north  of  25N  over  all  four 
meridians . 

Only  those  gridpoints  where  the  value  of  at  least  one  of  CL 
and  CL„  was  greater  than  0.1  were  included  in  those  populations.   The 
gridpoints  so  excluded  were  considered  climatologically  unrepresenta- 
tive for  the  purpose  of  cloud- tuning  the  albedo.   The  A (ALB)  data  was 
additionally  divided  by  season  giving  a  total  of  eight  separate  cases 
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further  defined  in  sequence  as  WIN-TR,  WIN-HL,  SPR-TR,  SPR-HL,  SUM-TR, 
SUM-HL,  AUT-TR  and  AUT-HL. 

The  Bimedical  set  of  programs  (Dixon,  1973)  was  utilized  to 
compute  the  root-mean-square  (RMS)  of  A (ALB)  values  corresponding  to 
a  wide  range  of  f-values  for  each  of  the  seasons  and  regions.   Values 
of  f  leading  to  the  minimum  RMS  of  A (ALB)  were  isolated  and  appear  in 
Table  VIII,  along  with  the  RMS  difference  of  A (ALB). 

TABLE  VIII.   Values  of  f  and  corresponding  minimum  RMS 
values  of  A (ALB) . 

Minimum  RMS 
f  of  A  (ALB) 

WIN-TR  .40  .069 

WIN-HL  .80  .073 

SPR-TR  .35  .067 

SPR-HL  .65  .082 

SUM-TR  .60  .079 

SUM-HL  .70  .109 

AUT-TR  .40  .052 

AUT-HL  .60  .073 

As  can  be  seen  from  Table  VIII,  the  listed  values  of  f  in 
winter  indicates  the  largest  difference  between  the  tropical  and  extra- 
tropical  cloud- types.   The  evenly  layered  cloud  formations  of  the 
northern  latitudes  in  mid-January  seems  more  closely  to  approximate 
those  considered  by  C.  D.  Rodgers  (1967) .   On  the  other  hand,  the 
smallest  difference  in  f-values  occurs  between  SUM-TR  and  SUM-HL 
indicating  a  less  significant  difference  in  cloud  types  over  the 
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complete  latitude  range  during  mid-July.   As  would  be  expected  the 
Spring  and  Autumn  f-ranges  fall  in  between  the  Winter  and  Summer  cases. 
A  full  listing  of  zonally-averaged  ALBMOD  values,  after  tuning  at  all 
gridpoints,  and  of  zonally-averaged  ALBRAS  values,  generated  from 
identical  gridpoints,  can  be  found  in  Table  IX. 

Upon  examination  of  Table  IX  it  is  apparent  that  the  tuned 
model-albedo  values  compare  closely  with  the  satellite  albedos  derived 
from  Raschke.   Table  IX  includes  the  annual  averages  and  the  Northern 
Hemisphere  cosine-weighted  means  as  well  as  the  mid-seasonal  averages. 

Further  inspection  of  the  weighted  means  shows  that  ALBMOD  is 
somewhat  smaller  than  ALBRAS  for  each  of  the  four  seasons  and  annual 
average.   However  most  of  this  underestimate  seems  to  occur  in  the 
HL-population  where  the  result  of  minimizing  A (ALB)  led  to  values  of  f 
substantially  less  than  unity  (Table  VIII) .   This  result  was  unexpected 
and  is  probably  due  to  a  mismatch  between  the  model  clouds  (Eq.  (2-2)) 
and  those  averaged  in  the  satellite  climatology  particularly  in  the 
higher  latitudes. 

The  conclusion  is  that  reflectance-tuning  to  satellite  albedo 
appears  valid  in  the  tropics  and  subtropics  where  cloud-formations 
on  a  given  day  tend  to  be  persistent.   However  tuning  by  this  method  in 
high  latitudes  is  less  conclusive  because  interdiurnal  cloud  variability 
is  much  greater. 

In  any  case,  all  further  computations  of  solar  dispositions  in 
this  study  will  utilize  the  appropriate  cloud-reflectance  tuning  factors 
as  listed  in  Table  VIII,  depending  upon  season  and  latitude. 
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V.   MERIDIONAL  CROSS-SECTIONAL  DEPICTION 
OF  THE  RADIATIVE-BALANCE  COMPUTATIONS 

A.  GENERAL 

The  general  design  of  this  section  is  to  utilize  all  of  the  compu- 
tational concepts  discussed  in  Sections  III  and  IV  in  the  computations 
for  a  single  time-step  in  the  radiative  heating  package  developed  for 
use  in  the  FNWC  prediction  model.   After  testing  ALBMOD  with  the  corres- 
ponding values  of  Raschke  et  al.  (1973)  ,  it  was  decided  that  only  the 
computations  by  the  appropriate  cloud-reflectance  tuning  factors  (Table 
VIII)  would  be  displayed  in  the  meridional  cross-sections  (Figs.  10... 
17)  and  in  the  comparisons  of  the  four-layer  versus  the  two-layer  cool- 
ing rates  summarized  in  Table  X. 

The  radiative  calculations  were  performed  at  each  gridpoint  of  the 
four  meridians  for  each  of  the  four  mid-seasonal  days  considered;  how- 
ever, only  the  winter  and  summer  cases,  as  illustrating  the  two  extreme 
situations,  are  presented  in  this  section. 

B.  GEOGRAPHICAL  REPRESENTATION  OF  THE  RADIATIVE-BALANCE  DISTRIBUTION 
The  FNWC  gridpoint  processed  analyses  for  0000GMT  were  used  at  the 

three  Pacific  cross-sections,  while  that  for  1200GMT  were  used  for  the 
single  Atlantic  meridian.  This  was  done  so  that  the  set  of  gridpoints 
was  considered  to  be  subject  to  the  actual  radiative- transfer  calcula- 
tions involved  for  these  specific  times.  Figure  9  depicts  in  symbolic 
language  the  key  to  the  computational  entries  in  Figs.  10,  11,... 16, 
17.  This  symbolic  list  presents  the  computations  made  at  each  radia- 
tive sounding  gridpoint  (I, J) having  data  in  the  form  of  Table  I.   The 
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computations  proceed  from  the  tropopause  (approximately  level  k=2)  to 
the  ocean  surface.   For  purposes  of  climatological  data  comparison, 
Figs.  10,  11,... 16,  17  were  developed  by  interpolating  gridpoint  results 
to  integral  multiples  of  5-degree  latitudinal  increments.   The  inter- 
polation routine  to  this  gridpoint  spacing  made  use  of  the  Lagrangian 
cubic  interpolation  scheme  (after  Spaeth,  1975) . 

bb)  -a,  "-""'?! "1"  +  o.  "-" (I-3) (I"4) 


1    (1-2)  (1-3)  (1-4)         x2     (2-1)  (2-3)  (2-4) 


(1-1)  (1-2)  (1-4)  (1-1)  (1-2)  (1-3) 

*•-}      n^n   /->_~>\   r~>—/\\  ~/ 


(5-1) 


•3     (3-1) (3-2) (3-4)  ~4     (4-1) (4-2) (4-3) 


Finally  for  ease  in  reconciling  the  magnitudes  of  all  radiative-trans- 
fer rates,  the  time-dependent  solar  disposition  rates  have  been 
averaged  to  24-hourly  rates. 

C.   EXPLANATION  OF  SYMBOLIC  TERMS 
1.   Cross-Section  at  Level  k=2 

The  discussion  of  all  insolation  parameters  discussed  previous- 
ly in  Section  IV  dealt  with  the  specific  time  of  day  that  corresponded 
to  the  hour  angle  h  for  the  instantaneous  time  t  under  consideration. 
The  incident  solar  insolation  dealt  with  is  then 


F(2)  =  S(— )  2  Cos  z  .  (5-2) 

r 
ra 


In  order  to  avoid  reference  to  specific  map  times  t,  the  instantaneous 
solar  hour-angles  were  h  =  35,  10,  55,  and  35  degrees,  respectively 
for  cross-sections  1,2,3  and  4  as  depicted  in  Figs.  10,  11,  12  and  13 
for  winter,  and  in  Figs.  14,  15,  16,  and  17  for  summer. 
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QAVE  represents  the  24-hour  average  of  F(2)  and  appears  as  the 
first  input  symbol  in  Fig.  9.   Its  value  is  considered  to  be  more  re- 
presentative climatologically  for  the  data  day  under  consideration  than 
F(2). 

QAVE  is  derived  by  the  formula 


2™=  "  F<2>  Hr!  <5-3' 


where 


Cos  z  =  [H  Sin<J)Sin6  +  Cosc})Cos6SinH]  /"IT  (5-4) 

H  =  ArcCos  [-Tan(}>TanS]  (5-5) 

Here  6  is  the  appropriate  solar  declination  angle  as  listed  in  Table  V, 
and  H  is  the  appropriate  hour  angle  at  local  sunset  at  latitude  (J).  The 
value  of  H  also  depends  upon  which  mid-seasonal  date  is  being  considered. 


Cos  z  in  Eq.  (5-4)  is  equal  to  the  24-hour  average  cosine  of  the  zenith 
angle,  Eq.  (4-2) .   The  24-hour  time  averaging  period  for  QAVE  gives  heat- 
ing results  consistent  in  magnitude  with  the  terrestrial  flux  divergences, 
which  change  only  slightly  with  the  time  of  day.   The  conversion  to  ex- 
pected daily  averaged  solar  disposition  quantities  is  compatible  with 
the  determination  of  a  hemispheric  radiative  balance  for  the  given  mid- 
seasonal  dates  (Sec.  IV.A.). 

Other  parameters  needed  for  level  k=2  are 


QREF  =  REF(t)  (      )  (5-6) 

Cos  z 


where 


REF(t)  =  F(2)  -  A24  -  A46  -  A68  -  A810 

-  (IA10  +  IS10)  (5-7) 
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REF(t)  is  the  instantaneous  solar  reflected  insolation  at  a  gridpoint 
and  QREF  is  its  24-hour  average,  assuming  that  the  instantaneous 
planetary  albedo  remains  constant  for  the  24-hour  period.   This  assump- 
tion requires  that  the  cloud  amounts  computed  at  the  indicated  synoptic 
times  are  representative  of  the  entire  day. 

The  same  principle  will  be  used  with  regard  to  all  other  solar 
parameters  in  the  conversion  from  time-dependent  values  at  solar  time  t 
to  24-hour  averaged  values.   Superior  bars  (   )  are  not  used  in  the  sym- 
bolism for  the  averaged  values  shown  in  the  cross-sections  key,  Fig.  9, 
but  are  implied  by  the  use  of  QAVE,  etc.  in  the  solar-disposition  terms. 
The  24-hour  average  tropopause  balance,  BALT,  is  computed  from 

* 

BALT  -  QAVE  -  (QREF  +  F   )  (5-8) 

for  the  level  k=2  at  the  indicated  latitude.   Net  terrestrial  fluxes, 

* 

such  as  F   ,  were  considered  to  be  constant  throughout  the  24-hour  period, 

a  valid  assumption  if  the  cloud  cover  remains  quasi-constant  for  the 
period. 

2.   Cross-Sections  in  Layers  (2,4) , (4 , 6) , (6, 8)  and  (8,10) 

The  following  definitions  apply  for  the  four  layers  identified 
in  the  present  section  heading  and  are  further  identified  in  Fig.  9. 
All  of  the  heat  transfers  shown  in  these  layers  are  assumed  to  be  of 
radiative  character  only.   The  daily-averaged  radiative  heating  (cooling) 
rate  in  layer  (2,4)  is  given  by 

BAL  24  -  024  -  F24  (5-9) 
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where 

Q24  =  daily-averaged  solar  absorption  in  layer  (2,4) 

and  is  defined  relative  to  A24 (t)  by  a  cosine  transformation  similar 
to  Eq.  (5-3) ,  and 

F24  =  terrestrial  cooling  rate  in  (2,4)  . 

Similarly  the  daily-averaged  radiative  heating  (cooling)  rates  in 
layers  (4,6),  (6,8)  and  (810)  are  given  by 

BAL46  =  Q46  -  F46  (5-10) 

BAL68  =  Q68  -  F68  (5-11) 

BAL810=  Q810-  F810  (5-12) 

respectively. 

3.   Cross-Section  at  Air-Sea  Interface  (k=10) 


The  radiative  balance  at  the  earth's  surface  (BALB)  is  as 
defined  in  the  following  equation: 


* 


BALB  =  QABG  -  F,_   .  (5-13) 

10 


QABG  is  the  24-hour  average  solar  insolation  absorbed  by  the  surface 
as  follows: 


QABG  =  QABG(t) (Cos  z/Cos  z)  .  (5-14) 

D.   MERIDIONAL  CROSS-SECTIONS  OF  THE  VERTICAL  RADIATION  BALANCE 

Figs.  10,  11,... 16,  17  as  previously  explained  represent  the 
single  time  step  of  heating  computations  for  each  of  the  four  meridians 
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used  in  this  study.   The  eight  figures  have  been  divided  into  (a)  tropi- 
cal results  and  (b)  mid-to-high-latitude  results  for  the  mid-winter  and 
mid-summer  cases  respectively.   While  the  results  depicted  in  these 
cross-sections  are  exhibited  as  representing  daily-averaged  values,  they 
are  actually  based  upon  radiation-computations  at  the  specific  map  times 
of  OOOOGMT  and  1200GMT  on  16  January  and  16  July  1974.   Therefore,  for 
these  results  to  be  meaningful  as  a  stepwise  part  of  the  heat  package 
subroutine  of  FNWC,  the  solar  radiative  absorption  and  reflectance  terms 
would  have  to  be  recoverable  as  a  function  of  GMT,  i.e., 


F(2,t)  =  QAVE*  (Cos  z/Cos  z)  (5-15) 


REF(t)  =  QREF*  (Cos  z/Cos  z)  (5-16) 

etc.   Thus  solar  disposition  terms  may  then  be  utilized  in  connection 
with  the  one-hour  stepwise  application  of  the  thermodynamic  equation  of 
the  set  of  primitive  equations  used  in  the  FNWC  prediction  process, 
assuming  the  "2/3-CL"  parameterization  of  Eq.  (2-2) . 

E.   COMPARISON  OF  THE  FOUR- LAYER  WITH  THE  TWO-LAYER  FLUX  DIVF:RGENCE  MODEL 

In  previous  studies  (A,  B,  C,  D)  which  considered  only  the  two-layer 
flux-divergence  model  [layers  (2,6)  and  (6,10)]  it  was  assumed  that  the 
radiative-cooling  rate  in  layer  (2,6),  BAL26,  was  uniformly  distributed 
over  layers  (2,4)  and  (4,6).   A  similar  remark  is  applicable  to  BAL610 
relative  to  the  sublayers  (6,8)  and  (8,10).   In  an  effort  to  determine 
the  vertical  resolution  of  the  atmospheric  flux-divergences  separately, 
the  four  values  of  BAL24,  BAL46,  BAL68  and  BAL810  have  been  averaged 
over  all  gridpoints  for  the  four  computational  dates. 
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The  computed  values  of  BAL24  and  BAL68  were  then  compared  with  one- 
half  of  BAL26  and  BAL610  respectively  to  determine  the  relative  per- 
centage of  heating  (cooling)  in  these  layers,  as  depicted  in  Table  X 
for  the  mid-seasonal  winter  and  summer  dates.   It  is  noted  that  the 
values  of  BAL46  and  BAL810  are  identical  in  magnitude  but  of  opposite 
algebraic  sign  to  BAL24  and  BAL68  respectively. 

The  results  of  these  four-layer  computations  (Table  X)  indicates 
that  in  layer  (2,4)  there  is  approximately  thirty  percent  less  cooling 
than  would  have  been  deduced  in  the  simpler  two-layer  model.   This 
implies  less  radiationally  induced  instability  in  the  upper  troposphere 
than  that  originally  estimated.   In  the  layer  (4,6)  there  is  compensat- 
ing increased  radiative  de-stabilization.   This  radiative  difference 
in  the  layer  (4,6)  appears  to  be  primarily  a  result  of  the  placement  of 
clouds  in  both  models.   For  instance,  the  upper  clouds  have  tops  at 
k=4,  and  IR  net  flux  acts  as  the  chief  source  of  increased  de-stabili- 
zation in  layer(4,6)  (and  of  decreased  de-stabilization  in  layer  (2,4)). 
Likewise  the  lower  clouds  have  tops  at  k=8;  so,  it  follows  that  IR  net 
flux  acts  as  the  chief  source  of  increased  de-stabilization  in  layer 
(8,10),  and  finally,  of  reduced  de-stabilizaticn  in  layer  (6,8). 

The  tuning  of  solar  reflectances  has  not.  introduced  a  significant 
radiative  influence  on  the  results  of  the  de-stabilization  differences 
just  outlined.   Rather  it  was  the  placement  of  clouds  relative  to 
layer  boundaries  that  affected  primarily  the  IR  net-flux  calculations. 
The  results  summarized  in  Table  X  are  then  to  be  considered  valid  if 
the  cloud-parameterization  model  approaches  climatological  reality. 
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VI.   THE  ZONAL  DISTRIBUTION  OF  RADIATIONAL 

BALANCE  TERNS  OF  THE  OCEAN-ATMOSPHERE  SYSTEM 

A.   GENERAL  INTRODUCTION:   ZONAL  CROSS-SECTIONS 

The  zonally  distributed  cross-sections  of  the  radiation  contributions 
over  the  ocean-troposphere  system  are  presented  in  Figs.  18,  19,  20  and 
21  for  winter,  spring,  summer  and  autumn  seasons,  respectively.   The 
cross-sections,  show  the  results  after  averaging  over  the  four  meridians 
on  the  mid-seasonal  dates  considered  in  this  study.   The  results  are  dis- 
played in  the  format  of  Fig.  9. 

In  obtaining  zonally  distributed  seasonal  means  of  radiative  heating 


rates,  denoted  0  ($)  ,    at  each  five-degree  multiple  of  latitude  (f)  in  the 
range  20S,...  65N,  values  of  each  radiative  parameter  listed  in  Fig.  9 
were  arithmetically  averaged  over  the  four  meridians.   At  <J)  =  65N  there 
was  only  one  contribution  to  the  zonal  average,  while  in  the  Southern  Hemi- 
sphere latitudes  20S,  15S  and  10S,  only  two  values  (on  X   =   125W  and  A  = 
35W)  of  each  parameter  contributed  to  the  means.   At  5S,  there  were  three 


meridional  sets  of  radiative  parameters  entering  into  Q((j))  .   Otherwise, 
there  were  four  seasonal  values  of  each  radiative  parameter  entering  into 


the  arithmetic  mean  Q(4))  at  each  latitude  (J)  of  the  cross-sections,  Figs. 
18 , . . . , 21 . 

Consequently,  near  the  northern  and  southern  boundaries  of  Figs.  18, 
...,21,  the  listed  seasonal  values  may  not  be  as  representative  as  in 
mid-latitudes.   Nevertheless,  the  general  equatorial-to-polar  trend  in 
the  radiative-change  terms  appears  reliable  in  both  hemispheres. 
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A  zonal  cross-section  depicting  the  mean  annual  radiative  distribu- 
tion has  been  constructed,  Fig.  22,  by  arithmetic  averaging  over  the 
four  mid-seasonal  sets  of  results  (Figs0  18,..., 21). 

It  has  also  been  convenient  in  earlier  sections  (cf.  Tables  III,  IV, 
VIII)  to  compute  "weighted  averages"  with  respect  to  latitudes  in  the 
Northern  Hemisphere.   Thus  a  definition  of  a  weighted-average  parameter 

which  takes  into  account  the  number  of  observations  k .  available  at  each 

l 

latitude  has  been  given  (Meyers,  1975)  as 


14 


i=l    4   \j=l   k.  / 


Cos  (J) . 


Q  =  Q  =    (6-1) 

*wt.  avg.       14 

ri.     ~  cos  (j). 

i=l    4       l 


Here  Q..  is  the  Q-value  on  meridian  j  at  latitude  d> .  and  k.  =  1.....4  is 
the  number  of  meridional  observations  available  for  the  arithmetic 


average  Q((f>.).   Note  that  i  =  1,...,14  corresponds  to  the  14  latitudes, 
<f> .  =  0 ,  . . . ,  6  5N . 

B.   ANNUAL  RADIATIVE  BALANCE  FOR  THE  EARTH-TROPOSPHERE  SYSTEM 

For  purposes  of  this  summary,  the  term  BALT  (at  k=2)  in  Fig.  9  has 
been  redesignated  for  simplicity  as  R  (0) ,  and  annual  values  of  R  have 
been  plotted  in  Fig.  23  as  a  function  of  latitude.   R  represents  the 
net  radiative- transfer  rate  across  the  tropopause,  that  is  at  the  top  of 
the  troposphere-ocean  system  considered  here.   Similarly,  the  net  flux 
at  the  surface  previously  denoted  BALB  in  the  key,  Fig.  9,  is 
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redesignated  here  following  Malkus  (1962)  as  R((J))  and  its  annual  distri- 
bution has  been  graphed  against  latitude  in  Fig.  23.  Finally,  the  rela- 
tionship between  R  and  R  in  any  column  is  given  by 


R     =   R  +  R  (6-2) 

t        a 


where  R   is  the  overall  net  cooling  rate  or  flux-divergence  of  the  tropo- 


spheric  column,  e.g., 


R   =  BAL24  +  BAL46  +  BAL68  +  BAL810  (6-3) 

a 


In  Fig.  23,  values  of  R  ((f))  have  also  been  plotted  against  latitude;  how- 

a 

ever,  a  simpler  method  of  obtaining  R  from  (6-2) ,  namely, 

a 

R   =  R^  -  R  (6-4) 

at 

has  been  employed.   The  resulting  zonal  annual  distributions  of  each  of 

the  three  parameters,  R  ,  R  and  R  ,  are  shown  as  functions  of  latitude 

t         a 

in  Fig.  23,  where  they  are  superimposed  against  similar  functions  drawn 
from  Sellers  (1965) „ 

The  comparison  which  follows  focuses  on  the  annual  distribution  in 
the  Northern  Hemisphere  only.   It  should  also  be  noted  that  Sellers' 
radiative  parameters  were  taken  from  climatology  at  the  top  of  the  mean 
zonal  atmosphere,  whereas  this  radiative  model  gives  corresponding  results 
for  level  k  =  2,  based  upon  tropospheric  soundings  over  an  oceanic  surface 
only. 

Throughout  the  latitude  range  (0-65N),  the  radiative  net  flux  R  at 
the  ocean  surface  is  substantially  greater  than  the  climatological  amount 
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Ra 
•  ^  RESELLERS) 

Ra  (SELLERS) 
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10  20  30 

LATITUDE 


50    60  65 N 


Figure  23.   Radiative  net  fluxes  at  the  tropopause  (R  )  at  the  ocean 
surface  (R) ;  and  net-flux  divergence  for  the  tropospheric  column  (R_) 
from  the  zonal  annual  results  (Fig.  22(a,b)).   Solid  lines  indicate 


values  of 


R  .  R  and  R  after  Sellers  (1965) 
t         a 
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of  Sellers.   This  is  attributable  to  the  reduced  model-albedo  values  of 

Sec.  IV. F.   This  in  turn  allowed  more  insolational-absorption  in  the 

ocean  thereby  increasing  R  over  all  latitudes  as  compared  to  Sellers' 

values  (Fig.  23) .   The  term  R-R   (which  equals  -R  )  may  be  regarded  as 

u  a 

an  atmospheric  "greenhouse  effect"  in  contributing  to  warming  of  the  sur- 
face both  for  the  model  and  for  Sellers'  climatology.   For  each  latitude, 
Fig.  23  indicates  that  the  "model  greenhouse"  is  in  close  agreement  with 
that  of  Sellers. 

The  R  -distributions  of  both  Sellers  and  of  the  present  model  are  in 
a 

good  agreement  across  the  entire  range  of  Fig.  23,  both  representing  cool- 
ing rates  in  the  troposphere  for  all  <J>.   A  minimum  value  of  R  at  latitude 
20N  appears  on  both  curves  and  is  presumably  due  to  a  local  maximum  of 
IR  net-flux  divergence  F*  -  F   *  (resulting  from  the  minimum  opaque 
cloud  cover)  in  subtropical  latitudes.   From  5 ON  to  65N,  an  increase  in 
the  mode]  values  of  R  apDears  in  Fig.  23  and  seems  to  indicate  a  reduced 
IR  net-flux  loss  associated  with  high  CL-values  (Table  III)  at  these 
northerly  latitudes.   Sellers'  results  indicate  the  opposite  trend,  i.e., 

increasingly  negative  R  -values  from  50-60N,  but  his  climatology  may  not 

a 

include  the  kind  of  detail  that  would  indicate  reasonable  oceanic  cloudi- 
ness over  these  latitudes. 

The  radiative  model  gives  R  -values  which  are  generally  larger,  both 
over  tropical  latitudes  as  well  as  over  the  northern  latitudes,  than  the 
R  -values  of  Sellers.   As  in  the  R((J))  comparison  this  is  a  result  primar- 
ily of  the  reduced  model-albedo.   In  the  mid-latitude  zone  25-40N  of 
Fig.  23,  close  agreement  occurs  between  the  two  R  -curves  presumably  due 
to  such  similar  c lima to logical  effects  as  for  example,  the  mean  polar 
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front.   The  positive-to-negative  R  crossover-point  occurs  very  near  to 
40N  for  both  R  -distributions  of  Fig.  23.   The  model  weighted-average  is 


■1  . 


R  ((f))  =  0.048  ly  min  '  if  averaging  is  considered  only  to  65N.   If  nega- 
tive R  annual  values  exist  in  latitudes  <p   >  65N,  an  ocean-troposphere 
radiative  balance  should  very  nearly  be  established  over  the  Northern 
Hemisphere.   However  data  for  the  model  were  not  analyzed  over  ice-covered 
regions  poleward  of  65N,  where  presumably  higher  surface  albedos  and 
cooler  drier  soundings  would  enable  R  to  assume  larger  negative  values. 

The  radiative-model  values  of  R  ((f)),  R((f>)  and  R  ((f))  are  in  agreement 

t  a 

with  present  albedo  and  IR  net-flux  observations  from  recent  satellite 
climatology.   It  should  be  noted  that  Sellers'  radiative  results  are 
based  on  climatology  of  London  (1957) ,  and  include  also  data  from  the 
land  areas  of  the  Northern  Hemisphere;  on  the  other  hand  the  radiative 
model  values  have  been  based  only  upon  soundings  over  the  ocean. 

C.   CROSS -SEASONAL  EFFECTS  IN  THE  NORTHERN  HEMISPHERE 
1.   On  the  Net  Flux  Across  the  Tropopause,  R 

The  seasonal  distribution  of  the  tropospheric  net  flux  is  pre- 
sented as  a  function  of  latitude  in  Table  XI (a)  where  the  listed  R  -values 
correspond  to  the  tropopause,  level  k  =  2.   One  may  discern  two  geographic 
regimes  of  annual  waves  in  Table  XI (a) ;  a  polar  zone  (25-65N)  and  a 
tropical-subtropical  zone  (0-25N) . 

The  cross-seasonal  variation  of  R  in  the  northern  regime  has 
the  following  property:  R  -values  are  a  minimum  in  winter  increasing  to 
a  maximum  in  summer  followed  by  a  decrease  in  autumn.   In  the  tropical 
zone  (0-10N) ,  R  -values  have  positive  peaks  in  spring  and  autumn.   The 
summer  minimum  of  R  in  0-10N  is  presumably  attributable  to  increased 
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cloudiness  associated  with  the  ITCZ  in  this  latitude  zone.   The  subtropi- 
cal region  15-25N  does  not  have  a  minimum  in  midsummer,  and  may  be  simply 

classified  as  intermediate  in  behavior  between  that  exhibited  by  R  in 

1      t 

the  polar  and  the  tropical  regimes.   These  geographic  cross-sectional 
effects  (Table  XI  (a) )  tend  to  confirm  the  general  validity  of  the  radia- 
tive model. 

2.   On  the  Sea-surface  Model  Balance,  R 

The  seasonal  distribution  of  the  net  radiative  flux  R  at  the 
earth's  surface  is  listed  as  a  function  of  latitude  in  Table  XI (b) .   Again 
a  polar  zone  (25-65N)  and  a  tropical-subtropical  zone  (0-25N)  are  defined 
for  discussing  cross-seasonal  effects  on  R. 

The  model  results  for  R  are  analogous  to  those  just  specified  for 
R  ,  namely: 

(i)   in  the  polar  latitudes  there  is  a  single  sine  wave  with 

annual  periodicity  and  maximum  in  midsummer; 
(ii)   in  tropical  latitudes  0-10N,  the  double  peaked  character 

with  maximum  in  spring  and  autumn  appears; 
(iii)   there  exists  an  intermediate  subtropical  zone  of  R-values 
which  no  longer  exhibits  a  minimum  in  summer. 

D.   TOP  OF  THE  ATMOSPHERE  COMPARISON  OF  MODEL  WITH  SATELLITE  DATA 
1.   Net  Radiative  Transfer  Rate  at  Top  of  the  Atmosphere 

The  three  parameters  that  describe  the  radiative  flux  at  level 
k  =  0  are  associated  in  the  following  way: 


R  =  p   -  F  •        (6-5) 

N    y-N  T 


Q  ,  the  daily  averaged  incoming  solar  flux,  can  be  related  to  QAVE  and 
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QREF  (Fig.  9)  by  the  relationship 


QN  =  |t||  QAVE(l-ALB)  =  |~|  QAVE  -  QREF  (6-6) 


where  QAVE  and  QREF  have  already  been  evaluated  according  to  Eqs.  (5-3) 
and  (5-6)  respectively  and  the  factor  2.00/1.92  has  the  effect  of  deriv- 
ing the  mean  daily  solar  insolation  at  level  k  =  0  before  ozone-oxygen 
absorption.   F  ,  the  IR  net  flux  to  space,  was  previously  discussed  in 
Sec.  III.D.  [Eq.  (3-11)],  and  R  ,  the  difference  between  Q   and  F  , 
represents  the  daily-averaged  radiative  net  flux  at  the  top  of  the 
atmosphere.   In  order  to  more  easily  distinguish  between  model-values 
of  R  and  the  same  parameter  from  satellite  climatology  of  Raschke  et  al. 
(1973) ,  R   is  further  identified  as  R  MOD  or  R  RAS,  respectively.   QAVE 
is  assumed  identical  both  for  the  model  and  for  the  treatment  of  the 
Raschke  data  in  this  study. 

2.   Zonally-Averaged  Computations  of  R 

With  the  use  of  Eqs.  (6-6)  and  (6-5),  values  of  Q  ,  F   and  R 

N    T      N 

for  each  of  the  four  mid-seasonal  dates  and  by  five-degree  latitude 
intervals  have  been  computed  both  for  the  model-parameters  and  the  satel- 
lite-observed (RAS)  parameters.   The  results  of  the  computations  both 
seasonally  and  for  the  annual  mean  case  are  presented  in  Table  XII. 

A  comparison  of  R  MOD  with  R..RA3  from  the  zonal-annual  results 
at  the  bottom  of  Table  XII  shows  that  the  crossover  from  positive-to- 
negative  values  of  R^  occurs  near  40N  for  both  model-  and  Raschke-results, 

A  statistical  accuracy  check  of  R  MOD  as  against  R  RAS  data  is 

N  N 

conducted  using  the  standard  deviations  of  the  difference  parameter 
(RMOD  -  R  RAS)  over  all  Northern  Hemisphere  latitudes.   A  summary  of 
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the  mean  and  standard  deviations  of  this  parameter  over  the  latitude 
range  of  these  five  cross-sections  (Table  XII)  follows: 

16  Jan.     16  Apr.     16  Jul.     16  Oct. 


RMOD-R  RAS      .004        .024        .019        .014 
Std.  Dev.      .028        .037        .053        .039 


Annual 

.015   ly 

.    -1 

mm 

.039    ly 

.    -1 

min 

The  zonally-averaged  mean  (R  MOD- R  RAS)  is  positive  for  each 
season  as  well  as  for  the  annual  distribution  (whose  mean  is  0.015  ly 
min   ) .   This  small  positive  difference  was  anticipated  from  considera- 
tion of  Sec.  IV. F.  (Table  IX)  where  the  middle- to-high  latitude  albedos 
were  tuned  slightly  too  small. 

The  overall  standard  deviation  in  (RMOD-R  RAS)  of  0.039  ly  min 

is  of  the  same  order  of  magnitude  as  that  of  Q  MOD-0  RAS.   This  compara- 

N     N 

tively  small  difference  in  the  case  of  model  versus  satellite  climatology 
was  a  considerable  improvement  over  (A,B,C,D)  where  untuned  albedos  were 
used. 


E.   STRATOSPHERIC  MODEL  RADIATIVE  BALANCE 

It  is  clear  that  the  flux-convergence  in  the  stratosphere  is  given 
by 


** 


P. 


where  R   is  given  by  Eg.  (6-5).   Note  that  R  is  identical  to  BALT  of 
N  t 

Fig.  9  which  reduces  to 


R^    =   OAVE  -  QREF  -  F  *  .  (6-7) 

t    "  2 
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Forming  the  difference,  R  -  R  ,  leads  to 


R  -  R  =  .04  OAVE  -  (F_  -  F  *)  .  (6-8) 

N    t  T     2 


The  stratospheric  flux-convergence,  R  -R  ,  is  displayed  in  zonal- 
annual  format  in  Table  XIII.   Based  on  the  Northern  Hemisphere  weighted 
results  (bottom  line) ,  it  appears  that  a  positive  net-flux  convergence 
exists  between  0-65N.   However  the  trend  in  R  -R  at  higher  latitudes 
is  towards  increasingly  negative  radiative  values  in  the  stratosphere 
poleward  of  65N  such  that  an  annual  mean  radiative  balance  may  be  in- 
ferred.  For  instance,  assuming  that  in  the  latitude  range  65-90N,  an 
average  zonal-annual  value  of  R  -R  =  -0.0502  ly  min  "  exists  at  77. 5N, 
then,  the  cosine-weighted  mean  is  -0.0126  ly  min    in  the  region  pole- 
ward of  65N.   Thus  the  summed  Northern  Hemisphere  weighted  mean  would  be 
zero  and  a  radiative  balance  would  exist  in  the  stratosphere. 

F.   ZONAL- ANNUAL  NORTHERN  HEMISPHERE  HEAT  BUDGET 
1.   Tropospheric  Radiation  Budget 

The  radiative  heating  rate  of  a  tropospheric  column  may  be  ex- 
pressed as  a  function  of  R  ((p)    by  the  right  side  of  Eq.  (6-9) 

a 

Q    +  S   =  R  +  (E  +  H„)  (6-9) 

xva    a    a         T 


S   is  the  storage  heating  rate  of  the  column  and  O    is  the  required 
a  va 

flux-convergence  of  sensible  and  latent  heat  energies  compatible  with 
the  heat  balance  at  latitude  (J).   R  is  the  tropospheric  radiative  net 
cooling  rate  and  is  depicted  in  Fig.  23.   (E  +  HL,)  includes  the  latent 
and  sensible  heat  parameters  respectively,  but  these  parameters  are 
not  part  of  this  study  and  will  be  considered  to  be  contained  in  Q 
of  Eq.  (6-9)  for  simplicity. 

Ill 


Lat. 

R 

N 

\ 

\-\ 

2  OS 

.1179 

.1030 

.0149 

15 

.1260 

.1085 

.0175 

10 

.1293 

.1074 

.0219 

5 

.1456 

.1248 

.0208 

0 

.1443 

.1215 

.0228 

5 

.1337 

.1097 

.0240 

10 

.1230 

.0969 

.0261 

15 

.1084 

.0826 

.0258 

20 

.0907 

.0644 

.0236 

25 

.0603 

.0367 

.0236 

30 

.0441 

.0293 

.0148 

35 

.0221 

.0113 

.0108 

40 

.0081 

.0052 

O0029 

45 

-.0137 

-.0107 

-.0030 

50 

-.0273 

-.0190 

-.0083 

55 

-.0426 

-.0168 

-.0258 

60 

-.0566 

-.0290 

-.0256 

65N 

-.0563 

-.0106 

-.0457 

Wt .  Avg . 

.0602 

.0476 

.0126 

TABLE  XIII.   C 

:omputed  radiative 

net  fluxes 

at  the  top  of 

atmosphere  (R  )  at  the  tropopause  (R  ) ;  and  the 
stratospheric  absorption  (R  -R  ) .   All  parameters 
are  listed  in  the  zonal-annual "format  in  ly  min~l 
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The  Northern  Hemisphere  cosine-weighted  R  turns  out  to  be 

R  =  -  0.1551  ly  min" 
a 

and  if  it  is  temporarily  attributable  to  mean  storage-cooling  S   of  the 

3. 

troposphere  alone,  the  mean  storage  rate  corresponds  to  a  temperature- 
change  rate  given  by 

ST         1440  (R~) 

(6t)  =  4#1  ~AP °C/day  (6-10) 

mb 


with  AP   =  800  mb  in  the  troposphere.   The  resultant  cooling  rate  (6-10) 

over  the  tropospheric  depth,  800  mb,  with  zero  lateral  flux  divergence, 

2 
is  approximately  -1.14°C  per  day  averaged  over  the  mean  cm  tropospheric 

column. 

Thus  the  general  circulation  of  the  atmosphere  would  have  to 

function  to  bring  about  the  atmospheric  heat  balance.   This  process 

would  require  the  proper  flux-convergence  Q    to  offset  the  annual  radia- 

2 
tive  loss  of  the  average  cm  column  of  the  troposphere. 

2 .   Zonal-annual  Heat  Budget  of  the  Ocean 

An  equation  similar  to  (6-9)  may  be  written  for  the  zonal-annual 

heat  budget  for  the  ocean  as  a  function  of  R  ((£>). 


Q    +  S   =  R  -  (E  +  Hp)  (6-11) 

vo       o  r 


Here,  S   is  the  storage  rate  of  the  ocean  water  mass  and  0   is  the 
o  3  *vo 

required  mean  oceanic  heat  flux  divergence  for  a  balance  at  latitude  (J). 
The  primary  ocean  mass  heating  parameter  is  the  net  radiative  heating 
flux  at  the  ocean  surface,  R  (Fig.  23),  and  (E  +  Hp)  is  for  convenience 
considered  part  of  Q 
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The  Northern  Hemisphere  annual  weighted  mean  value  of  R  is 


R  =  0.2027  ly  min 


which  corresponds  (Table  XI)  to  a  mean  radiative  heating  rate  in  the 

water-mass  column.   It  should  be  noted  that  this  mean  heating  rate  is 

noticeably  larger  than  R  and  of  opposite  sign  to  the  corresponding 

a 

tropospheric  radiative  cooling  effect  (R  =  -0.1551  ly  min   ) „   How- 

a 

ever,  if  the  three-dimensional  oceanic  flux  divergence  of  both 
sensible  and  latent  heat  is  considered,  the  oceanic  heat  balance  should 
result. 
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VII.   CONCLUSIONS 


The  major  change  in  this  radiative  model  from  the  model  evolving 
from  (A,B,C,D)  was  the  introduction  of  tuned  cloud  reflectances.   This 
innovation  brought  the  resultant  model-albedos  into  agreement  with 
satellite  climatology  for  comparable  data  dates. 

Gridpoint  comparisons  at  the  top  of  the  atmosphere  of  model  net  flux 

F   and  of  the  net  incoming  model  insolation  Q  were  made  with  similar 
T  N 

parameters  from  Raschkes '  satellite  climatology  (1973).   Closer  agree- 
ment was  obtained  between  Q  MOD  and  Q  RAS  than  in  any  of  the  preceding 

N         N 

studies  (A,B,C,D) .   Similarly,  close  agreement  in  a  least-squares 
sense  was  achieved  here  between  F  MOD  and  F  RAS  primarily  because  of 
an  improved  formulation  for  the  F  -computation  than  was  used  in 
earlier  studies.   Finally  the  ocean- atmosphere  system  net  flux 

N  "  CN  "   T 

has  only  a  small  least-square  error  when  model  and  satellite  results 
are  compared  with  properly  stratified  gridpoint  and  time  cross-sections. 

In  this  study,  mid-latitude  cloud  reflectances  were  also  tuned  to 
obtain  closer  agreement  between  the  global  albedo  and  corresponding 
values  from  satellite  climatology.  The  result  was  to  systematically  in- 
crease Q„  relative  to  Raschke's  values.   This  result  was  undoubtedly 
N 

due  to  a  mismatch  between  the  model  clouds  and  those  involved  in  the 
climatological  average.  The  conclusion  reached  is  that  if  tuning  of 
cloud  reflectances  of  this  nature  is  to  be  attempted  in  mid-latitudes 
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using  FNWC  soundings,  more  accurate  cloud-parameterization  formulas 
should  be  established  against  observed  cloud-cover  amounts  for  identical 
experiment-periods  similar  to  GATE.   If  then  further  cloud-reflectance 
tuning  is  necessary,  it  could  proceed  from  a  more  certain  knowledge  of 
the  "ground  truth"  provided  by  these  experiments. 
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APPENDIX  A.      COMPUTER  PROGRAM 

//FGRT.SYSIN    CD    * 

COMMON    67 ( 11) ,DA7A(11,5) ,CL(2) 

COMMON/ ABS/    F A2 , AL  PH  A G , C Z , S ECZ , F S , FA  D J 

COMMON/TRIG/    S7  ,  CT,SD,CD,  ALAT 

CCMMCN/VP/    Q(ll) ,QS(11) 

CCMMCN/ARM/    RM 

CCMMON/RASKC/RASVAL 

CGMMC N/RASK/    AL BRAS (93) ,F2RAS(93) ,  ISO 

CCMMGN/FACT/    K,IS£A 
502    FQFMAT(4X,2(F8.4,2X) ) 
500    FGPMAT     ( 2X , F6. 1 , IX , F7. 2 , 1 X , 5( F 9. 5 , IX) , 6X , 21 2) 

DC    68    IS2A=1,4 

IFdSEA.EQ.il    DECa-21. 0791 7/57 ,29578 

IF( TSEA.-Q.2)    DEC  =       8.48333/    57.29578 

IF  dSEA.EQ.3 )    DcC=       21.5000/57.29578 

IFdSEA.EQ.4)    DEC=         -8,225/57.2  9576 

IFdSEA.EQ.l)    RM*.  9837 15**2 

IFdSEA.EQ.2  )     RM=1. 00333**2 

IF(  IStA.EQ.3)    PM=1. 01644**2 

IFUSEA.EQ.4)    RM*.  997165**2 
C 

C       7FE    FGLLOWING    READ    STATEMENT    APPLIES    TO    THE    SATELLITE 
C       CLIMATOLOGY    DATA    G^    RASCHKB       ET    AL.     AS    USED    IN    THE    Sol 
C       ROUTINE    REFT. 
C 

DC  9  ISC=1,93 

READ(5,  502)  AL BR* S ( I  SO ) ,F2RAS ( I  SO ) 
9  CONTINUE 

CC  6  6  ISO =1,93 

CO  10  1=1,11 

FEAD(17,500,END=99)     (DATA  (I  ,K),K  =  l,2),0d),CSd), 
1 (CATA( I  ,  KK) ,KK=3,5I ,LAB1,LAB2 
1C    CONTINUE 

I F  dS  E  A.  EQ. 3.  AND.  I  SG.EQ.  25)     LAB  I  =  24 

CALL    ANGLE(CZ,SSCZ,ISO,LABlfDEC> 

IFdSEA.EQ. LAND. ALAT. LE. 25.)    RASVAL    =.4 

IFdSEA. EQ. LAND. ALAT. GT. 25.)     RASVAL    =.6 

IF(ISEA.EQ.2.AND.ALAT.LE.25.)    RASVAL    =.35 

IFdSEA.EQ. 2. AND. ALAT. GT. 25.)    RASVAL    =.65 

IFdSEA. EQ. 3. AND. ALAT. LE. 25.)     RASVAL    =.6 

IF(ISEA.EQ.3.AND.ALAT.GT.25.)    RASVAL    =.7 

IFdSEA.EQ. 4. AND. ALAT. LE. 25. )     RASVAL    =.4 

IFdSEA. EQ. 4.  AND, ALAT. G7. 25.  )     RASVAL    =.6 

C  i L  L    EBB 

ALPFAG=(.54*(.7-CZ) )+.06 

IF( ALPHAG.L7..C6 )     ALDFAG=.06 

FS=.651*CZ*1.92/RM 

FACJ=(2.00*CZ)/RM 

FA2=( .349*1. 92*CZ)/RM 

DO    20    1=1,4 

CL(1)=0.0 

CL(2)=0.G 

IFU.EQ.2.3R.I.EQ.4)    CL(1)=1.000 

IF(I.GE.3)    CL(2)=1.000 

W7=(1.G-DATA(6,3) ) *( 1 .O-DATA ( 2, 2 ) ) 

IFd.EQ.2)    WT=DATA(6,3)*d.0-DATA<2,  3)  ) 

IFd.EQ.3)    WT=DATA(2,3)*d.0-DATA<6,3)  ) 

I F ( I . EQ . 4 )     WT=  DA T  A ( 2 , 3  )  -DAT A ( 6,3) 

CALL    A6SuRB(CL(l ) ,CL(2),WT) 
20    CONTINUE 
6  6    CONTINUE 
66    CCNTINlE 
99    STOP 

ENC 
G 

SLERGUTINE    ABSORB ( CI , C2,WT ) 

CCMMON/ABS/    FA2, AL PHaG, C Z , S ECZ , FS, FADJ 

COMMON/ R AS KC/RASVAL 

IF(C1.GT..5.0R.C2.GT..5)    GO    TO    20 

ALPHAR=(.085+< . 25074* ALuG 10 ( SEC Z ) ) ) 

IF(ALPHAR.GT.1.000»     ALPHAR= 1.00  0 
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A=(  <1.-ALPHAR)*(1.-ALPHAG)  )/(  1 . - < ALPHAR* ALP  FAG )  ) 
SOL=A*FS 

GC  TC  2  1 


SLBROUTINE  TD< ANS, tND,Nl,N2 ,SEC,C) 

CCMMON  ET( 11) ,DATA(  11, 5),CL<21 

U  =  DATA(M,4)-  DATA(N2,4) 

A=.303 

ANS=C*<  (U*SEC  )**A  ) 

IF( IND.GT.IO. )  ANS=1.G-ANS 

RETURN 

ENC 

SLBROUTINE  BBB 

CCMMCN  BT(  11  )  ,  DAT A (11, 5) ,CL<2) 
SIGMA =1.1 70403 
E=l.E-7 
CO  ID  1=1,11 
A=CATA( 1,2 ) +2 73. 16 
ET(I )=( A**4)*SIGMA*E 
10  CGMINUE 

CATA(2,4)=(DATA( 1 , 4 > +OATA ( 3 ,4 ) ) /2 .00 
DATA(2,5)=(DATA( 1 , 5 ) +DATA ( 3 , 5) }/2.C0 
RETURN 
ENC 

SLBROUTINE  ANGLE (CZ,SECZ,IS,L2,D) 

CCMMQN/TRIG/  ST, CT , SO, CO, AL AT 

P=57. 29578 

h=35.0/R 

IFCIS.LE.25)  H=55.0/R 

IF (IS .GE.26.AND. IS.LE.50)  H=1C.0/R 

A=573.752 

S=2.0 

I  F  (  I S  . G  E  .  26.  A N'D  .  IS.LE.50)  S  =  1.0 

AK=L2+1 

IFUS.GE.68)  AK=53-L2 

B=(S*U32.-AIO**2)   ) 

S7=(A-6)/( A+B) 

CT=(1.-(ST**2) )**.5 

SC=SIN(C) 

CC=COS(D) 

CF=CGS( H) 

CZ=(ST*SD)+(CT*CD*CH* 

SECZ=1.00/CZ 

TANT=ST/CT 

ALATMATAN(TANT) )*57. 29573 

RETURN 

ENC 

SLBROUTINE    AB00(WT,SQL) 

CGMMGN /  A  b-S /  F  A2 , ALPHA  G , C  Z , $  EC  Z , F  S , FAD J 

C£LL    TD(U24,0,9, 7,SECZ,.271  ) 

CALL    TC(U26,0,9,  5,  SSCZ,.2.71) 

C*LL    TC(U23,0,9,3,SECZ,.27I  ) 

CALL    TC(U21G,0,9,1,SECZ,.271) 

CALL    TC(U1C, 12,9, 1,SECZ, .271) 

A24=FA2*U24 
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A46=FA 
A68=FA 
A610=F 

TRANA= 
AI10=T 

PEFA=F 
CALL  R 
RETURN 
ENC 


2*(U26-U24) 

2*(U28-U26 ) 

A2*(U210-U28) 

FA2*U10 

RANA*( l.-ALPHAG) 

A2-A24-A46-A6  8-A81C-AI 10 

EFT( 1, A24,A46, A63,A81C, M 


10,TRA«MA,REFA,S0l_tWT  ) 


SLB 
CCM 

ccr 

Rl  = 

Al  = 
C*L 
C^L 
CAL 
CAL 
FAD 
F4U 
F2L 
A24 
A  46 
F6C 
F8C 
F10 
F10 
F&'u 
FfcU 
F6C 
F8C 
F10 
A6£ 
A61 
All 
TRA 
RfcF 
CAL 
RET 
ENC 


2,9,7,SECZ,.271) 
2t5»3»l»6667» .271) 


ROUTINE    AB10< W7,SQL)  , 

M0N/ABS/FA2tALPHAG,CZ,SECZ,  FS,FADJ 

MON/PASKC/RASVAL 

•46*RA$VAL 

.20 

L    TC 

L    TO 

L    TO 

L    TD 

=  FA2 

=R1* 

=  F4U 

=  FA2 

=  F4D 

=  F4D 

=  F6C 

D=F6 

U=Fi 

=  F1C 

=  F10 

D=F6 

C=F6 

DC-F 

=  F60 

0-F8 

0=(F 

NA  =  F 

A=FA 

L    RE 

URN 


024,1, 

D68?l, 

0  61.0,12,5,1,  I.0666 7,  .271) 

0310 

D24 

D 

024 

4D  +  F41 

1 

1.-R1 

068 

TD610 

*ALPH; 

TD810 

TD610 

Rl 


12  ,3,1, 1.6667, .271) 

U-F2U 
-Al  ) 

I 

IAG 


0*10  66 

0 

c 


0*70610 

8D+F8U-F6U+F6DD-F8DD 
-F10D+F10U-FSU+F8DD-F10DD 

00+{ l.-ALPHAG) )+F10DD 

0D+F10DD 

-A24-A46-A68-A81G-AI 10 

T(2, A24, A46, A6  8, A81Q, A  I 1C ,TR ANA , REF A, SOL , WT ) 


SLBROUTIN 

CCNMCN/AB 

CCMMON/hA 

R2=.50*RA 

A2=.30 

CALL    TD( 


E    A301(WT,S0L) 

S/FA2,4LPH*G,CZ,SECZ,FS,FA0J 

SKC/RASVAL 

SVAL 


CALL 
CALL 
CALL 
C  ALL 

C^LL 


TC( 
TC< 
T0< 
TC< 
T0( 


024,12,9,7,SFCZi.271 

026 
028 
091 


,12,9,5,SE 


Li 


F4D^FA2*T 
F6C=FA2-T 
F6C=FA2*T 
F£U=F8C*R 

FfcL=F8U*~ 
F4L=F8U*T 
F2l=F8U*T 
A24=FA2-F 
A4fc=F4C-F 
A68=F6C-F 
F9C=F8D*( 
F1Q0=F9C* 
F1CU=F1C0 
F9L=F1GU* 
F9C0=F9L* 
F1GDC=F9D 
A  8 10= (F 80 
TRANA=F10 
AI10=(F10 


048 
024 
026 
D23 
2 

068 

048 

0  28 

40+ 

60  + 

30^ 

1.- 

T09 

*AL 

709 

R2 

D*T 

*A2 

OF 

D*( 


c71) 
,i^,9,3fS5CZ, .271) 
0,12  ,2,1, 1.6667, .271) 
,12, 5,3,3tCZ, .271) 
,i2,7,353ECZ, .271 ) 


F4U-F2U 

F6U-F4U 

F8U-F6U 

32-A2) 

10 

PHAG 

10 

0910 

)+F9  0-F10D+F10U-F9U+F9DD- 

1000 

l.-ALPHAG) )+F10DD 


•FIODO 
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PEFA=FA2-424-A46-A68-A810-AI10 

CALL    RcFT<3,  A24 , A4  6 , 46  8, A3  10,  A  I  10,TRAN4,  REFA,S0L,W7) 
RETURN 


SL5 
CCN 
CCM 
Rl  = 
R2  = 
Al  = 
A2  = 
CiL 
C£L 
CAL 
F4C 
F4U 
F2U 
A24 
A46 
F6D 
F8C 
F8L 
F6L 
F6C 
F8D 
A68 
F9C 
FIG 
FIO 
F9L 
F9D 
FIO 
A81 
All 
T*A 
REF 
CAL 
RET 
ENC 


ROUT 
MON  / 
MC  N  / 

•  46* 

.50* 

•  20 
.30 

L  TC 
L  TC 
L  TC 
=  FA2 
=F1* 
=  F4U 
=  FA2 
=  F40 
=  F4D 
=  F6C 
=  F8C 
=  F8U 
D=F6 
C  =  F6 
=  F6C 
=  F8D 
C=F9 
U=F1 
=  F10 
U=F9 
DC=F 
0=(( 
0=(F 
Na=F 
A=FA 
L  RE 
URN 


INE    A611(WT,SCL) 

A8S/FA2,ALPHAG,CZ,SECZ,FS,FADJ 

PASKC/RASVAL 

RASVAL 

P.ASVAL 


(TO 
(TO 
(TO 
^.=  TD 
F4D 
*TD 
-F4 
*A1 
*  (  I 
*TD 
*R2 
*TD 
U*R 
DD^ 
-F8 
*(  1 
D*T 
00* 
U*T 
U*R 

9  DO 
F80 
100 

1 0  D 
2-4 
FT< 


24,12,9.  7,SECZt.27D 
66,12, 5,3,1.6667,  .271  ) 
910, 12,2,1, l.o667, .271  ) 

2h 

24 
D+F4U-F2U 

•-R1-A1) 

68 


68 
1 

TD 

0  + 

C9 

AL 

D9 

2 

^7 

+F 

*( 

+  F 

24 

4, 


68 

F8U-F6U+F6DC-F8DD 
R2-A2)+F8DD*( 1.-A2) 
10 

PHAG 
10 

D910 

8DD)*A2 J+F90-F1CD+F 1GU-F9U+F VDD-F 10DD 

l.-ALPHAG) )+F10DD 

100D 

-A46-A68-A810-AI1Q 

A24, 446, A68,A810, AI 10,TRAN4, REFA,SOl  ,WT) 


StePOUTINE    REFKNSUB,  A24,  A46,  A68,  A810,  A  I  10,TRANA, 
lREFA,SOL,wT) 


10 


CO 

cc 
cc 
cc 
cc 
cc 

IF 
CA 

CA 
CA 
CA 
CA 
CS 
CP 
CA 
CA 
CA 
CA 
CA 
CS 
CR 
CA 

IF 
AE 
TA 
TA 
TA 


MMON 

MMON/ 

NMCN/ 

MiMCN/ 

M.MON/ 

PMCN/ 

MKCN/ 

(NSU8 

24=0. 

46=0. 

63=0. 

610  =  0 

1 10  =  0 

110=0 

EF  =  0. 

24  =  Ca 

46  =  C4 

68=CA 

810  =  C 

I10  =  C 

I10  =  C 

EF  =  CR 

LL    IR 

(NSUB 

G=(CA 

ND=SC 

NT=ST 

N  =  -TA 


ET 

4B 

TR 

AR 

RA 

PA 

FA 

•  G 

0 

0 

\J 

.0 

.0 

.0 

Q 

24 
46 
68 
A8 
A  I 
SI 
EF 
(N 
.L 
II 
/C 
/C 
NO 


(11  ),  DATAUli  5)fCL<2) 

S/    F42, ALPHaG,CZ,5ECZ, FS.FAOJ 

IG/    ST,CT,SD.CD,ALAT 

M/    RM 

SK/    ALBRAS(93} ,F2RAS(93)  ,ISG 

SKC/    RASVAL 

CT/    K.ISE-A 

T.l)    GO    TO    10 


+  A2  4* 

■♦■446* 

+  A68* 
10  +  48 
10  +  AI 
10  +  SO 
+  ((FS 
SUB,  A 
T.4) 
0  +C  S I 
D 
I 
*TANT 


WT 

WT 
WT 
10 
10 

L* 

- o 
10 
Re 
10 


W ' 


w  i 

OL  )**T)+RcFA*WT 

, A3,A6,A4,A2,F10,F8,F6,F4,F2,rtT,FT) 

TURN 

)/FADJ 
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H  = 
SI 

BA 
F2 
QA 

AF 
QP 
EA 
C2 
AF 
BA 
Q4 
AF 
BA 
Q6 
AF 
BA 
08 
AF 
BA 
CA 
BA 
AF 
CA 
C2 
AF 
EA 
CA 
Q6 
AF 
BA 
CI 
DI 
DI 
DI 
AL 
AL 
F2 
RN 
RK 
CI 
RE 
EN 

SU 

lhT 

IF 
OF 
El 
E£ 
E6 
54 
E2 
ET 
CA 
CA 
CA 
CA 
CA 
CA 
El 
E8 
E6 
E4 
E2 
ET 
IF 
El 
E& 
E6 
E4 


ARC  OS 
NH  =  SI 
RCCS  = 
00=1. 
VE=F2 
2=-F2 
EF=(C 
LT  =  C 
4=(CA 
24=-( 
L24  =  G 
6=<CA 
46  =  -( 
L46=C 
E  =  (CA 
68=-( 
L  68=  Ci 
10=(C 
610=- 
L310  = 
EG=  ( 
LB=  G 
1G=-F 
26=CA 
§=(•5 
26=  — . 
L26  =  C 

6iO  =  C 

*  n-  t 

j.  **  i  • 

610=- 
L610= 
F24  =  B 
F46  =  6 
F  6  8  =  B 
F810= 
E=CRE 
BDIF  = 
CIF  =  F 
RAS  =  ( 
MOD=< 
FRN=R 
TURN 
C 


(TAN) 
N(H) 

( <S7*SD*H)+(CD*CT*SINH)  )/ 2. 1415  9 
92/RV 
00*8ARCOS 


R 
A 
2 
F 

F 

4 

6 

F 

6 
j 

{ 

Q 

A 
1 

2 

5 

2 
f, 

H 

5 

G 

A 
A 

B 
F 
AL  t;- 


5F 

VE 

4* 

£  — 

4  + 

6* 

4- 

6  + 

8* 

6- 

8  + 

81 

F8 

81 

BG 

BG 

0 

4  + 

CA 

*( 

6  + 

68 

*C 

5* 

61 

L2 

L4 

L6 

AL 

/F 

LB 


*QAVE)/( FADJ*.9600) 
-QREF-F2 

QAVE) /(FADJ-,9600) 

F4) 

AF24 

QAVE )/<FADJ*.9600) 

F6) 

AF4  6 

3AVE)/(FADJ*.96C0) 

F8  ) 

AF68 

C*0AV6)/(FA0J*.9600) 

-F10) 

0 + A  F  8 1  0 

*QAVE) /.9cOQ 

-FIG 


CA 
26 
F2 
AF 
+C 
A  6 
(F 
C  + 
4- 
6- 
8- 
31 
AD 


AV 
AV 

MO 


2R 

C:': 

D- 


46 

*  Q  -    V   : 

-F6) 

26 

A  8  10 

LC*Q 

6-F1 

AF61 

BAL2 

3AL2 

BAL6 

O-BA 

J 

LERA 

ASU 

1.C4 

1.04 

RMR.A 


)/(FADJ*.96G0) 


AVE)/(FADJ*.9600) 
0) 

0 

o 

6 

10 

L610 

S(ISO) 

SO) 

1667*(1.-ALBRASUS0))  J-F2RASUS 

lbo7*(l.-AL8) )-FT 

S 


BROUTINE    IR(NSUB,FlQ,F8,F6,F4,F2,ElC>,E8,r6,E4,E2, 
tET) 

(NSUB.GT.l »    GO    TO    10 
F  =  1.0 
C  =  0  .  c 
=c.o 

=  0.0 

=  0.0 

=  0.0 

=  0.0 

LL    FF10(F10,0FF) 

LL    FF8(F8,0FF) 

LL  FF6(F6,0FFi 

LL  FF4(F4,0FF) 

LL     FF2(F2tOFF) 

LL    FFTGP(FT,OFF) 

0=E10+WT*F10 

=E8+WT*F3 

=E6+WT*F6 

=E4+WT*F4 

=E2+WT*F2 

=ET+WT*FT 

(NSUB.lT.4)    return 

C=E10/1440. 

=E8/1440. 

=E6/1440. 

=E4/1440. 
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E2=E2/1440. 
ET=ET/1440. 

RETURN 
ENC 


1C 


SL 

IF 
CA 
CA 
CA 
CA 
C/i 

c  a 
C/s 
El 
B2 
Al 

ei 

62 
E3 
64 
A  2 
Bl 
A3 
Fl 
1  (C 
RE 
EN 

SL 
CC 

I  F 
CA 
CA 
CA 
CA 
CA 
C* 
61 
62 
E3 
64 
85 
Al 
A2 
A3 
F8 
1*A 
RE 
EN 

SL 

CC 
IF 
CA 
CA 
CA 
CA 
CA 
CA 
Bl 
6  2 
B3 
34 
85 
Al 
A2 


10    Ft 
1*A 

RE 


ERG 

cc.v 

(CL 

LL 

LL 

LL 

LL 

LL 

LL 

LL 

=  (6 

=  (6 

=  (C 

=  BT 

=  (B 

=  (B 

MS 

=  (i 

=  (6 
C=( 
L(2 
TLR 
C 


UTI 

>rjN 
(l  } 

EWC 
EWC 
EWC 
EWC 
EWC 
EWC 
BCU 
7(  3 
T{  1 
BTC 
(10 
T(  9 
TC7 
T(  5 
(6T 
.-( 
T(  1 
(1. 
)=M 
N 


3 


NE    FFlCt 

BT( 
.GT, 
(  EW8 
C  EW6 
(  EW4 
(EW2 
(EW1 
(  EW9 
NDCW 
)-BT 


11) 
.5. 
tCF 
iOF 
,OF 
»0F 
,GF 
tOF 
(AVE 


T 

F» 

~  » 

r  t 
,0 
)* 


(5) 
)-BT(3)  ) 
1)-BT(5)  )- 

VE 

(  10 


10 
AT 

3, 

5, 
7. 
9. 
10 

2, 
FF 
(E 

EW 


Fl 

1,5) ,CL(21 
).GT..5)     Gi 


.*WA 
)-BT 
J-BT 
)-3T 

(5)  - 
•  5*E 

)-67 
-CL( 
3) 


.OF 

AC  1 

LC2 

1) 

1) 

1) 

1) 
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